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ABSTRACT 

We test analytic predictions from different models of magnetospheric accretion, 
which invoke disk-locking, using stellar and accretion parameters derived from mod- 
els of low resolution optical spectra of 36 T Tauri stars (TTSs) in NGC 2264 (age~3 
Myrs). Little evidence is found for models that assume purely dipolar field geometries; 
howeve r, strong suppo rt is found in the data for a modified version of the X-wind 



model (|Shu et al.lll994l ) which allows for non-dipolar field geometries. The trapped flux 
concept in the X-wind model is key to making the analytic predictions which appear 
supported in the data. By extension, our analysis provides support for the outflows 
predicted by the X-wind as these also originate in the trapped flux region. In addition, 
we find no support in the data for accretion powered stellar winds from young stars. 
By comparing the analysis presented here of NGC 2264 with a similar analysis of stars 
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in Taurus (age~l-2 Myr), we find evidence that the equilibrium interaction between 
the magnetic field and accretion disk in TTS systems evolves as the stars grow older, 
perhaps as the result of evolution of the stellar magnetic field geometry. We compare 
the accretion rates we derive with accretion rates based on U-band excess, finding good 
agreement. In addition, we use our accretion parameters to determine the relationship 
between accretion and H/3 luminosity, again finding good agreement with previously 
published results; however, we also find that care must be used when applying this 
relationship due to strong chromospheric emission in young stars which can lead to 
erroneous results in some cases. 

Subject headings: accretion, stars: low-mass, stars: magnetic field, stars: pre-main 
sequence, stars: variables: T Tauri 



INTRODUCTION 



T Tauri stars (TTSs), first classified and named bv lJovl (|1945l ). are low mass, pre-main sequence 
(PMS) objects that are roughly grouped into two classifications: 1. classical T Tauri stars (CTTSs), 
which show evidence of a circumstellar disk and mass accretion onto the central star in the form of 



exces s emission in the X-ray, UV, optical, and infrared (|Bertoutlll989l : lArgiroffi. Maggio. &: Pere 



20071 ). and 2. weak-line T Tauri stars (WTTSs), or naked T Tauri stars (NTTSs), which are also 
PMS stars bu t do not show evidence for significant mass accretion and do not have inner disks 
(jWalterl 1 1 98 71 ) . WTTSs are thought to be the evolutionary product of CTTSs that have ceased 
accreting material from their disks and are continuing their gravitational contraction to the main 
sequence, though there is significant overlap in HR diagrams. Understanding how CTTSs interact 
with and ultimately disperse their disks is vital to our knowledge of how these systems evolve, 
supposedly into the WTTS phase, and eventually become sun-like stars, some of which are likely 
surrounded by planetary systems similar to our own. The goal of this study is to test predictions 
of magnetospheric accretion theories (discussed below) in the ~3 Myr old open cluster NGC 2264. 
The slightly greater age of NGC 2264 compared to well-studied regions like Orion and the Taurus- 
Auriga star forming region, both with an age ~l-2 Myr, allows an exploration of magnetospheric 
accretion in an older population of TTSs which permits an investigation of how magnetoshperic 
accretion may evolve over time. 

It is now well established that the excess emis sion observed in the spectra of CTTSs is due to 

The 



the accretion of gas from a circumstellar disk (e.g. iBertout et al.lll988l : lHartigan et al 



1991 



origin al proposition for how the excess emission forms was put forth by iLvnden-Bell &: Pringle 
) who imagined the interaction between star and disk occuring in a boundary layer: the 
fast-rotating disk dissipates its energy in a narrow region which is in contact with the surface 
of the more slowly rotating central star, resulting in the excess emission observed at blue and UV 
wavelengths. While the boundary layer par adigm is able to explain the hot continuum emission (e.g. 
Bertout et al.1 1 19881 : iBasri fe Bertoutl Il989l ) , it fails to account for the large velocity-shifts in the 
line profiles of CT TSs and the large equivalent widths of Ha emission that are commonly observed 
(jHartmannl 1 19981 ) . Further evidence against the boundary layer model includes observed inner- 
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disk holes aroun d some CTTSs (jMever et al.lll997l ). which nevertheless shows signs of acc retion 



Hart mannlll998r). Chromosphe ric models explaining the observed excesses were proposed by [Cram 



19791 ) and ICalvet et al.l (jl984T ). and these were successful in reproducing some of the spectral 



chracteristics of CTTSs. However, like the boundary layer paradigm, these models fail to account 
for velocity shifts and widths in spectral lines and could not reproduce near-IR (disk ) excesses, as 



well as becoming increasingly unfeasible when used to describe highly veiled spectra (ICalvet et al 



19841 ). Magnetospheric accretion, in which gas from the disk is loaded onto stellar magnetic field 
lines and impacts the surface of the st ar at near free-fall v elocities, provides the best explanation 
for the observed properties of CTTSs (jBouvier et al.1 120071 ). 



Support for magnetospheric accretion as the dominant accretion proce ss on CTTSs is strong 



Magn etospheric accretion onto astrophysical objects was first investigated bv lGhosh. Lamb. &: Pcthick 



1977) who modeled pulsating neutron stars accreting mass from a binary companion. 



Uchida fc Shibata 



1984 ) first suggested a similar accretion mechanism for TTSs. The observational link between 



young stellar objects (YSOs) and the role magnetic fields play in their evolution was first estab- 
lished based on images of strong, collim ated bipolar outflows and m ore highly collimated jets, and 
the energy required to power them (see lAppenzeller k, Mundtlll989l ). 



The first investigations of mass loss from TTSs, b ased on ob servations of energetic winds, were 
performed by IVarsavskyi ([1960) , iHerbigl ()196ll ) , and iKuhil (119641 ) by assuming that the flows were 
driven from the surface of the star. Some of the first evidence for the circumstellar origin of outflows 
in CTTSs, and as a result the conclusion t hat these outflows mig ht be intrinsically linked to mass 
accretion onto the star, was discovered by Edwards et al.1 (|1989i ). They suggested that the winds 
and outflows were direct results of accretion and thus could not b e purely stellar in origin. We now 
know that bipolar outflows are ubiquitous phenemona in YSOs (jBallv. Reipurth. Davis 1 120071 ). 
These outflows require the existence of strong, hourglass-shaped magnetic fiel d lines to transport 
and collimate material from the disk or sta r into the observed bipolar flows (jPudritz Norman 
1983 



1986; 



Shu et al 



1988 



Camenzindlll990i ) ; winds driven purel y by therma l or radiation pressure 



cannot account for the large outflow velocities that are observed (|Ladalll985l ). Magnetic field lines 
can also act as a collima ting agent, bounding the outflows at large distances from the central 
star (e.g. IShu et al.lll988l ). Magnetospheric accr e tion theories can also account for th e observed 
variability of the UV excess (IBertout et al.l Il988l : lAlencar. Johns-Krull. fe Basril l200ll ) : accretion 
columns terminating at the stellar surface co-rotate with the star, thus producing variability on the 
same timescale as the stellar rotation period. Velocity shifts of emission and absorption lines are 
natu rally explained by the ~250 km/s velocities of the infalling material (e.g. ICalvet &: Gullbring 
19981 ). 



One of the most important questions concerning the evolution of TTSs is the problem of how 
thes e objects are able to shed angular momentum and rotate with velocities well b elow break-up 
(e.g. IVogel fc Kuhilll98ll : iHerbst et aLlbond : Lamm et allhoO-d : iMakidon et aDhoO-d ). The magne- 
tospheric accretion model provides a potential answer to this question due to the role of the stellar 
magnetic field in transferring angular momentum from the star to the disk: if the magnetic field 
couples with a sufficiently ionized disk and acts as a braking torque on the star, it will rotate more 



slowly than if no braking mechanism were present (|Konigllll99ll ; IShu et al.lll994l ). This interaction 
essentially locks the star to the disk (i.e. disk-locking) and prevents the star from rotating at 
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break-up velocity. 



This scenario can accoun t for many of the ob served rotation rates of PMS stars ([Edwards et al 



1993 



Kearns &: Herbst 



1998 



Lamm et al.l |2004| . hereafter L04) . Strong evidence for disk-locked 



stella r rotation has bee n found in the Orion Nebula Cluster (ONC) ( Choi Sz Herbstfl996 : Herbst et al 
20oi) and NGC 2264 dLamm et allboOjj hereaft er L05) , thou gh contradictory re sults have b 



)ccn 



found by the studies of lStassun et al.l (jl999T ) and iRebullI (|200ll ) for the ONC, and Makidon et al 



(|2004l ) for NGC 2264. Evidence su pporting magnetic disk- locking manifests itself in the form of 



a bimodal period distribution (e.g . 



Herbst et al 



Attridge fc Herbstlll992l ) and the detection of disk signatures 



(Edward s et al.lll993l ; iRebull et al.ll2006l ). a nd hence circumstellar disks, around long-period stars 



(2002 



1)^ based on the work of Choi Sz Herbst (1996), point out that the shorter pe- 



riod peak in the bimodal distribution is simply a binning artifact (i.e. the rotational evolution of 
stars that have ceased to be regulated by their disks should span a range of sh orter periods) while 
the longer period peak (~8days) is a direct resul t of magnetic disk- l ocking (lAttridge &: Herbst 



1992; Herbst et al 



2001 



200 2). To date, studies bv lRebull et all pOOff t . iDahm. Slesnick. fc White 



20121 ) , and ICieza fc Baliberl (|2007l ) provide the strongest eviden ce for correlations between slowly 
rotating stars and the presence of a circumstellar disk, though the lDahm. Slesnick. h White! (|2012l ) 
study finds high co nfidence levels (>99%) for only the M-dwarfs in their sample. The work of 
Herbst et al.l (J2002J) also shows strong evidence that longer period stars have a higher incidence of 
circumstellar disks. These correlations support the hypothesis of star-disk interactions being the 
main culprit in removing angular momentum from TTS systems. 



L05l also report a bimodal period distribution for a large sample o f stars in NGC 2264, similar 
to what has been reported for the ONC. However, iMakidon et al.l (120 04) find no evidence for 



a bimodal period distribution in NGC 2264. In fact, IMakidon et al.l (|2004l ) find no significant 
difference between the peri od distributions of the Orion region and NGC 22 64, w hich differ in 
age by a factor of ~2 (|L05l ). This result directly conflicts with the results of lL05l who find that 
the period distribution of angular momentum conservation for NGC 2264 is consistent with st ellar 
contraction from the period di stribution of the ONC based on fully convective PMS models. IL05 
explain that the conclusions of I Makidon et al.l (120041 ) are a result of the latter st udy c omparing the 
period distribution of NGC 2264 with a larger, inhomogeneous region in Orion. lL05l . on the other 
hand, only compare the period distribution of NGC 2264 with the younger homogeneous region of 
the ONC. Studies of both clusters find that the rotation period distributions for higher mass stars 
(M* > 0.25 Mq) and lower mass stars (M* < 0.25 Mq) peak at different locations, with the lower 
mass stars peaking at a shorter period than those with greater mass (jHerbst et al.l boOll . bood : 
Makidon et aDbood ; lLamm et aDbood ). 



Despite some conflicting results, observational evidence seems to point to disk-star interactions 
as being the primary candidate for explaining TTS rotation rates. Further support for magnetic 
disk-locking can be found in the me asurements of magnetic fie l ds on TTSs. Th o ugh t h e sample size 
of mea sured fields is relatively sma lk Johns-Krull et al. (j 19991 ). Guenther et al. (1999), Johns-Krull 
(|2007l ). and lYang Johns-Krulll (|201ll ) measure relatively uniform values of ~l-3 kG for surface 



magnetic fields on TTSs. The relative constancy of their magnetic fields, combined with the 
ubiquitous presence of circumstellar disks, points to interactions between the stellar magnetic field 
and disks as being a prime candidate for angular momentum regulation in CTTSs. However, 
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Johns-Krulll (|2007l ) points out that the magnetic fields of TTSs may not be strong enough to 
enfor ce disk-locking if a dipolar fiel d geometry is not assumed at the stellar surface. On the other 
hand. I Johns-Krull &: Gaffordl (|2002l ). hereafter JG02, investigated correlations predicted by several 
different theories of magnetospheric accretion (see below) and found support for models that assume 
a non-dipolar surface field geometry and only weak support for models assuming purely dipolar 
fields. 



I n this s tudy, we extend the analysis performed bv lJG02l to a larger sample of PMS stars in NGC 
2264. |JG02| examined several analy tic relat i onshi p s predicted by four different magnetospheric ac - 
cretion theories, specifically th ose oflKonigll (119911). ICameron Sz Campbelll (j 19931 ) . IShu et al.l (| 19941 ). 
and a modified version of the lOstriker &: Shul (|1995l ) model (hereafter OS95), using multiple sets 
of observations of CTTSs in the Taurus-Auriga molecular cloud complex. The theories and their 
predictions will be discussed in Section 2. Mode ling the accretion c olumns producing the excess 
emission in CTTSs as hot slabs of hydrogen (e.g. IValenti et al.lll993l ) enables us to derive reliable 
estimates of stellar and accretion parameters for our sample in NGC 2264. We then test the same 
correlations predicted by the aforementioned theories for this set of stars. NGC 2264 is an ideal can- 
didate for this study due to the availability of rotation periods and spectral type determinations for 
a large number of stars. As mentioned abo ve, N GC 2264 is slightly older than the Taurus-Auriga 
region and the ONC (3 Myrs vs. 1-2 Myrs; |L05| ) and so should provide a good testing ground for 
the models using a more evolved PMS stellar population. Exploring young clusters at different 
ages allows us to see when the equilibrium magnetspheric accretion relationships and disk-locking 
break down. Section 3 describes our observations and reductions and the classification of the stellar 
sample into accreting and negligibly accreting stars. In Section 4 we discuss our models and how 
they are used to derive stellar and accretion parameters (see Appendix A for a detailed description 
of our fitting procedure and model assumptions) . Section 5 includes the application of the accretion 
theory predictions to our sample and tests of the resulting correlations. A discussion of the results 
is given in Section 6 and our conclusions are presented in Section 7. 



2. ACCRETION THEORY PREDICTIONS 



Johns— Krull et al.l (119991 ) exa mined the magne tspheric accretion theories of iKonigll (|199ll ). 
Cameron fc Campbell ( 19931 ). and Shu et al. ( 1994 ) deriving equations for the stellar magnetic 
field (B*) in terms of i?*, P ro t, M, and M* where R* is the stellar radius, M* is the stellar mass, 
Prat is the rotation period of the star, and M is the mass accretion rate onto the star. One of the 
main assumptions underlying these accretion theories, and, as a result, the derived relationships, 
is that the central star is magnetically locked to its disk so that the star rotates at the Keplerian 



velocity of the disk t runcation radius (jKonigl 



1991 



Cameron Sz Campbell 



1993 



Shu et al. 1994 



Ostriker fc Shulll995l ). That is, the star and the inner truncation point in the disk are co-rotating. 
As a result of torques acting through the magnetic field, the central star is spun down and is not 
able to rotate at break-up velocities, whereas high rotation rates would normally be expected if 
the star is allowed to contract and accrete without angular momentum being removed from the 
system. By looking at the correlations predicted between stellar and accretion parameters by these 
magnetospheric accretion theories we can test the validity of these models as applied to PMS stars. 
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Of the models examined, the IKonigll (|199ll ) and ICameron fc Campbelll (|1993l ) theories are 
jlv mass accretion models and do not consider mass loss via a disk or stellar wind. The 



Shu et al. 


( 


1994 


) and 


OS9EI 1 



the central star with the launching of a magnetocentrifugal disk wind from the disk truncation 
point. The inclusion of wind-launching has important implications for the amount of matter that 
is ultimately accreted onto the star and also for t he amount of m agn etic flux that par ticipates in 



the accretion flow. The relationships predicted by IKonigll (|199ll ) and IShu et al.l (119941 ) depend in 
the same way on the stellar and accretion parameters: 



RlB div oc M^M^Pl/ 6 



rot ' 



(1) 



The lCameron &: Campbelll (|1993l ) prediction deviates only slightly from eq. (1): 



R*Bdi p oc 



^2/3^23/40^29/24^ 



(2) 



In eqs. (1) and (2), B^ip is the equatorial dipolar field stre ngth at the stellar surface. Plots 
created with eqs. (1) and (2) produce nearly identical results ( Johns-Krull &: Gafford 20021) . Fo r 
this reason we only d i splay plots generated using eq. (1) when considering the IKonigll (|199ll ). 
Cameron fc Campbelll (|1993l ). and IShu et al.1 (| 19941 ) models. 



Standard magnetospheric accretion models all assume that the stellar magnetic field is purely 
dipolar. Observations of magnetically sensitive abs orption lines in TTSs, however, have shown that 



most TTSs have complex surface fie l d geometries dJohns-Krul 



2004 



Daou. Johns-Krull. fc Valenti 



2006 



Donati et al 



200 



et al 



2011a 



1999 



Valenti fc Johns-Krull 



JG02 



show that if the 

dipolar field requirement is relaxed, but that the amount of magnetic flux participating in the 
accretion flow is conserved and mapped back to the stellar surface from the truncation p oint (this 
is equal to one-t hird o f the total magnetic flux trapped at the truncation radius; see IOS95I for 
details), then the lOS95l model can be modified to predict the following relationship: 



ncc-Bacc OC 



(3) 



where f acc is the filling factor of the accretion columns on the stellar surface, i.e. the percentage 
of the surface covered by accretion flows. In eq. (3) B acc is no longer the dipolar field strength 
but rather the strength of the field participating in the accretion flow at the surface of the star. 
The magnetic flux participating in th e accretion flow is equ al to one-third of the total magne tic 
flux trapped at the tr uncatio n point (jOstriker &: Shul I1995T). Equation (3) (eq. 7 from |JG02j) is 
a specific prediction of |JG02r s modified analysis of the IOS95I theory of magnetospheric accretion. 
More recently, investigations into the consequences of multipole magnetic field configurations on 
the observational properties of accretions flows have been undertaken using numerical models (e.g. 



1 OS95 is actu ally a det a iling of the accretion flow and the role of flux pinching at the disk trunc ation point 
predicted by the IShu et al.1 (|l994f ) model and not a new theory unto itself. Thus the IShu et al.l (|l994) and QS95 
references can be assumed to refer to the same model. 
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Mohantv Shull2008l ; LLong. Romanova. <fe Lovelace] 120081 ; I Adams fc Gregory||2012l ). providing fur- 
ther evidence that surfa ce magnetic fields on TTSs are probably not purely dipolar. However, 
Adams fc Gregory! (|2012l ) find that contributions from higher order magnetic field components are 
negligible at typical disk truncation radii in TTS systems. Thus the assumption that the dipole 
component of the magnetic field dominates when determining the truncation radius remains valid. 

Measurements of mean surface magnetic fields on CTTSs have shown the fie l d strengths to be 



relat i vely constant from s tar to star (IJohns-Krull et al. 



2004 : |johns-KrullH2007l : lYang fc Johns-Krullll201lh 



1999; 



Guen ther et al.lll999i ; IValenti Sz Johns-Krull 



Donati et al.l (|2011bl ) have suggested, how- 



ever, that the strength of the dipole and octopole components of the magnetic field on TTSs can 
vary strongly from star to star depending upon the existence and extent of a star's ra diative core: 
fully convective stars host strong, stable dipolar fields (e.g. iDonati et al .1 120081 . l2010l ) while stars 
with small radiative cores (0.0 < Mm rg /M*< 0.4] tend to have weaker di polar fields and mor e dom - 
inant octopolar field components ( Donati et al. 2011al lbl). For example, Donati et al. ( 2008 . 201C) 
measure 1.2 and 2-3 kG dipolar magnetic field components on the fully convective CTTSs BP Tau 
and A A Tau, respect ively; for the partially convective CTTSs V2129 Oph and V4046 Sgr, however, 
Donati et al.l (|2011al ) measure 900 and 50-100 G dipolar magnetic field components, respectively, 
for each star. Although the number of stars for which individual components of the magnetic field 
have been measured is relatively small, the Donati et al. results suggest that these field components 
vary from star to star. As mentioned above, the strength of the dipolar component has important 
implications for the location of the disk truncation radius. If the dipolar field is weak enough, the 
accretion disk may pinch the field closer to t he star allowing disk m aterial to accrete along more 



complex regions of the star's magnetosphere (IRomanova et al.ll201ll ). Such a scenario would most 



likely not be described by the equlibrium relationships of eqs. (1) and (3), due in large part to the 
fact that disk-locking can no longer be assumed if the truncation radius is much smaller than the 
co-rotation radius. 

Equations (1) and (3) are directly testable if estimates for the necessary parameters can be 
obtained. We note that f acc is not the same as the filling factor of the magnetic field on the 
surface of the star, which is taken to be 1.0 in the models being tested. Estimates for / acc , which 
are significantly les s than one, can be obtained by modeling the accretion flow as a hot slab of 



emitting hydrogen (IValenti et al.1 Il993l : lHartigan et al.Ml99ll ). This procedure will be discussed 



in §4. Due to the lack of magnetic field measurements for stars in NGC 2264, we are not able 
to substitute measured values of the surface magnetic field into eqs. (1) and (3). We instead 
proceed by first making the simple assumption that the dipolar component of the magnetic field 
(eq. 1) and the magnetic field participating in the accretion flow (eq. 3) are constant from star 
to star. Because we are concerned only with how well the two sides of eqs. (1) and (3) are 
correlated, assuming a constant field strength allows us to ignore B^ip and B acc in the analysis. 
While this assumption is justified based on measurements of mean surface magnetic fields on TTSs 
(see above), recent work by Donati et al. showing variations in the dipolar field component from 
one TTS to another potentially renders the constant field assumption inadequate, at least for stars 
which have developed radiative cores. In an attempt to roughly mimic this behavior, we use the 
internal structure calculations from the Siess et al. (j200ol ) PMS evolutionary models to determine 
the extent of the radiative cores in our sample. We then assign an equatorial dipolar magnetic field 
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strength based on the (limited) literature results in order to test the full forms of eqs. (1) and (3). 
This procedure and the results of both methods (assuming a constant field vs. assigning individual 
field values) are presented in §5. 



3. OBSERVATIONS AND DATA REDUCTION 



Our sample consists of 36 pre-main sequence stars of spectral type K and M in the open 
cluster NGC 2264 (a=06:41, 5=+09 35), and 14 main sequence dwarfs ranging from spectral type 
Kl to M 4 wh ich serve as spectral templates. The sample is a subset of o bjects from the rotation 
study of |L0J. The stellar identification numbers are those used by |L04j. In the end, only eight 
of the observed templates are actually used in our models. The NGC 2264 targets were chosen 
based on the availability of a determined spectral type and rotation period, and their confirmed 
PMS evolutionary nature based on approximate Ha equivalent widths. All of the stars in our 
sample have an Ha-index (a measu re o f the strength of Ha e missi on) greater than .1, which 
indicates the stars are likely CTTSs (|L05l ). We briefly discuss the |L04| classification of our sample 
in §3.3. The rotation periods were determined based on photometric variability studied by L04; the 
sp ectral types were de termined to within 1 subtype using optical and near infrared spectroscopy 
by lRebull et al.1 (|2002l ). All of the stars in our sample are brighter than I<16.0. 



3.1. Observations 

A summary of the observing runs is given in Table 1. We obtained low resolution (R~600) 
optical (3000 A- 6000 A) spectrophotometry of each of the stars listed in Tables 2 and 3. The 
wavelength resolution of the McDonald observations is ~7.2 A; for the Ki tt P eak observ ations 



AA~3.2 A. Spectrophotometric flux standards taken from lHamuv et al.1 ( 19921 ) and lOkel ( 19901 ) were 
observed each night at various airmasses in order to generate atmospheric extinction curves and 
flux conversion factors. Typical exposure times for the McDonald sample were ~1000 seconds with 
a signal to noise ratio of ~90-100 at 4400 A, approximately the center of the modeled wavelength 
region and free of any significant emission lines. Exposure times for the Kitt Peak sample were 
typically ~1200-1800 seconds with a signal to noise ratio of ~40-50 at 4400 A. Signal-to -noise 
ratios for fainter Kitt Peak stars (I> 15) are closer to ~20. Typically the S/N ratio of the data 
decreases towards bluer wavelengths as the detector becomes less sensitive to these photons. 

On January 29, 2004 light cirrus clouds were present during the observations of HD10476, 
L5575, 5967, 6039, and 5394. Exposures of L5967 and 6039 were retaken on Jan. 30 in order to 
obtain more accurate flux measurements. L5355, 1316, 3636, and 6024 were also observed through 
light cirrus on March 1, 2005. Calculated flux values for HD10476, L5575, 5355, 1316, 3636, and 
6024 should therefore be taken as lower boundaries to the actual stellar flux. This uncertainty is 
propagated through our model calculations and the final derived parameters for these stars are 
more uncertain than those derived for the stars with observations taken on photometric nights. 
These values should be viewed with caution. 
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3.2. Data Reduction 

All of the data were reduced using custom IDL routines. Bias and flat lamp images were 
taken at the beginning of each night. For the McDonald sample, argon comparison lamp spectra 
were taken at regular periods during the night in order to account for wavelength shifts along 
the detector; for the Kitt Peak sample, iron-argon lamp spectra were used. Cubic and quadratic 
dispersion solutions were calculated for each lamp spectrum using ~15 lines from 3300A-5500A. 
The lamp spectrum taken closest in time to the stellar exposure was used as the wavelength solution 
for that particular star. Median bias images were constructed and subtracted from each science 
image. In order to remove slit illumination effects by the flat lamp from the McDonald sample, 
polynomial fits were calculated perpendicular to the dispersion direction and then divided through 
the specified row in the flat image. These normalized flats were then combined to a create a master 
median flat for each night. The normalized fiat was divided into each science image in order to 
remove pixel-to-pixel variations on the CCD. Bad pixels and cosmic rays were removed manually 
by averaging over the adjacent two pixels in the dispersion direction. Each science image was then 
corrected for atmospheric extinction and multiplied by the flux conversion factor calculated for that 
night. The CCD read noise and statistical noise were included in the uncertainty calculations of 
the flux for each star. 



3.3. Grouping the Data 

Typical spectra and model fits are shown in Figure 1 (see Appendix A for the full sample), 
where the flux is given in units of 10~ 16 erg cm~ 2 s _1 A -1 . The observed spectrum is shown by 
the solid black line; the fit to the data is overplotted with a red line. The excess emission from the 
slab (dashed blue line) and the underlying scaled template (solid green line) are also plotted. Our 
sample is separated into three categories: strongly accreting stars, weak accretors, and negligible 
accretors. The negligible accretors (NA) were identified due to the lack of any excess emission 
required to produce good fits to the data, i.e. only a scaled and reddened template star was needed 
to match the observed spectrum. The inclusion of a slab in these cases yielded accretion rates of 
< 1CT 10 Mq yr _1 , almost two orders of magnitude lower than the average rate for the strongly 
accreting (SA) stars, and changes in the \ 2 goodness-of-fit measure of <10%. Model fits of the 
NAs are not shown in Figure 1. The weak accretors (WA) display negligible Balmer jumps but 
required a slab in order to reproduce the observed flux blueward of the jump. Chi-sqaured values 
for the WAs changed by >10% when a slab was included in the model, indicating that low levels 
of mass accretion are occurring and producing a small amount of excess emission. Figure 1 clearly 
shows that the slab emission is a larger percentage of the total emission for the SA stars than for 
the WA stars. Most SA stars show clear Balmer jumps and exhibit strong Balmer emission lines, 
as well as strong Ca II lines and some Fe n emission. Fits to the strongly accreting sample changed 
significantly (Ax 2 >50%) when a slab was included in the model. The average accretion rate for 
the SA stars is 10 times larger than the average rate for the WA stars. 

The distinction between strongly accreting stars and weak accretors is defined here quantita- 
tively by the ratio of the mass accretion rate to the stellar mass: if M/M* > 1 (in units of 1CT 8 
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yr -1 ) the star is categorized as a strongly accreting star; if M/M* < 1, the star is grouped as a 
weak accretor. While this is an arbitrary distinction, the groups are useful when examining the 
relationships predicted by eqs. (1) and (3). In addition, this criterion appears to separate the 
sample nicely into rough evolutionary regimes, with the majority of the SA stars having ages < 5 
Myr and the majority of weakly accreting stars having ages > 5 Myr. Age determinations for the 
sample are discussed in §4.4. The negligible accretors display a range of ages and these stars are 
most likely WTTSs whose disks have largely dissipated. By determining which of the stars belong 
to the NA category we are able to eliminate objects that might yield unrealistic slab and accretion 
parameters due to their very low or non-existent accretion rates. The disk-locking theories being 
tested here most likely no longer apply to these stars. 

The sample was chosen based on the classification of almost our entire sample as CTTSs 
according to a definition based on i?-band and Ha magnitudes. This definition is similar to the 
ori ginal classifica tion of the stars based on Ha equivalent widths, with CTTSs having W(Ha) > 5 
A ( Bertoutl ll989). However, it is obvious from the spectra of the negligibly accreting stars (Figure 
1, panel c) that these stars, and a small number of the weakly accreting stars, do not show typical 
char acteristics of CTTSs, e.g. strong Balmer emission lines and large veilings at blue wavelengths. 



L05l state that this classification definition selects stars that "most likely show the properties of 



accreting PMS stars" but that it does not provide strong proof of their inclusion in one group or 
another, due to the potentially large equivalent widths of Ha from chromospheric emission. Thus, 
we find that our SA and WA stars are probably CTTSs, while our NA stars are most likely WTTSs. 



4. MODELING THE EXCESS 



The excess emission from CTTSs is believed to result from an accretion shock at the stellar 
photosphere: matter accreting along magnetic field lines slams into the stellar surface at near 
free-fall velocities, i s strongly shocked a nd heated, and emits in the X-ray, UV, and optical (e.g. 
Kenvon et al.lll994l ; iKastner et al.ll2002l ). This excess emission is believed to be produced over a 
range of temperatures and densities along the shock column. The picture becomes more complicated 
if one considers the heating of the underlying photosphere, spot effects, and the interplay between 
the pre-shock and post-shock region. While simple 1- D shock codes have been utilized as a more 



realis tic method to produce the excess spectrum (e.g. ICalvet &: Gullbringl Il998l : iMuzerolle et al 



2003), uncertainties in the temperature structure of the accretion flows and how they are heated 
remain unresolved issues for more advanced shock models. For our purposes, however, all that is 
needed is a reliable estimate of the excess accretion emission, distinct from the stellar photospheric 
emission. 

Good estimates of the excess emission can be obtained by treating the accretion flow as a 
single-temperature accretion column, or slab, in LTE originating at the photosphere of the star 



(Valenti et al 



1993 



Gullbring et al.lll998l ). Esti mates derived from slab m odels show good agree- 



(Calvet h Cullbrine 


1998 


) and slab models con- 


young objects (e.g. 


Herczeg & Hillenbrand 


2008; 



Herczeg. Cruz. Hillenbrand! l2009l ). Although this is almost certainly not the physical scenario 
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producing the excess emission, order of magnitude estimates of the emission are good enough to 
the test relationships predicted by the accretion theories being investigated here. In any case, sys- 
tematic uncertainties in the mass accretion rate may dominate compared to differences produced 
by differences in the model treatment of the accretion flow itself, at least as far as slab and 1— D 
shock models are concerned. We discuss some of these systematic uncertainties in §4.3. 



4.1. Slab Models 



We model our sample of PMS stars by combining a scaled main sequence template star spec- 
trum with the scaled slab emission which represents the emission from the accretion column. The 
parameters governing the slab emission are number density, n, temperature, T, the length, / through 
the slab, and the filling factor, f acc , of the total stellar surface area covered by the accretion columns. 
The emi tting slab in the mod eled region (3400 A- 5500 A) is taken to be composed of pure hy- 
drogen. lHartigan et al.l (|199ll ) found that metals contribute <5% of the electrons in the slab for 



typical temperatures and densities in our models. In addition, lHartigan et al.l (|1989l ) show that 
emission line flux from all elements constitutes <5% of the total excess flux. For our purposes 
here we are primarily interested in estimating the accretion luminosity and resulting accretion rate, 
and uncertainties associated with the exclusion of elemental emission lines are acceptable since the 
uncertainties in the accretion rate are likely dominated by larger uncertainties, such as degeneracies 
in the fits themselves, geometric factors, and uncertainties in the distance and stellar radii. Thus, 
emission lines from other elements are not included in the model calculations. We also fix the 
turbulent velocity of the slab gas at 150 km/s in order to roughly match the line widths seen in 
our sample. The turbulent velocity actually serves as a substitute for the wavelength resolution 
of the instrument and does not represent a physical property of the slab. Letting the turbulent 
velocity vary in our fits does not result in significant improvements to the fit. If the line width is 
not well matched by Vr=150 km/s then Vr is manually adjusted until a good fit to the line width 
is achieved. 

The slab emission is calculated by solving the Boltzmann and Saha equations for the given 
temperature and density and then evaluating the optical depth and source function, in this case 
the Planck function, at each wavelength. The maximum number of levels included in the hydrogen 
atom (typically ~80) is determined by the temperature and density of the slab. The total opacity 
of the slab is calculated using the sum of Hi, Hn, H~, and electron scattering opacities. The most 
important feature of the slab is the Balmer jump, which occurs at ~3650 A and we define here 
to be (/b-/r)//r, the flux just b lueward of the jum p minus the flux just redward of the jump 
normalized to the red continuum (jValenti et al.lll993l ). Our low resolution spectra do not resolve 
the tightly packed higher Balmer lines, which merge into a pseudo-continuum around 3700 A. 
Thus, the Balmer lines blueward of H<5 are ignored in our fitting procedure. H7 and H/3 are also 
ignored w hen fitting the data du e to the possibility of winds significantly contributing to the line 
strengths (jAlencar &: Basrill2000i ). For this reason we have chosen to use only H<5 as a constraint on 
the Balmer line strengths. This generally produces good fits to the higher Balmer lines but poorer 
fits to H7 and H/3. However, we are mainly concerned with fitting the excess continuum and thus 
are not overly concerned with emission models that produce the best fits to the lines. To calculate 
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the emergent spectrum, we solve the radiative transfer equation through the slab. 

Our model procedure attempts to match the observed spectrum, particularly the Balmer jump. 
It is valuable to discuss how variations in the slab parameters affect the shape of the slab emission. 
The size of the Balmer jump is governed by the temperature of the gas and the optical depth 
through the slab. Increasing the temperature for a fixed density results in a weaker jump due to 
a faster increase in the populations of levels n>2 compared to increases in the n=2 population. In 
other words, the optical depth blueward of the jump decreases relative to the optical depth redward 
of the jump and the jump is diminished. The optical depth through the slab is controlled by n, 
I, and, to a lesser extent, T. Typical values for n and I that we find in our models are 10 14 cm -3 
and 10 7 cm, respectively. For these values and a slab temperature of 9000 K, the slab continuum 
is optically thin (r < 1) and the Balmer lines are very strong, i.e. r>>l in the lines. As the 
density is increased, the slab becomes optically thick and the line emission is weakened relative 
to the increasing continuum emission. For densities above 10 15 cm -3 , the slab becomes optically 
thick and is dominated by H~ continuum emission. Increasing I will also increase r. Because our 
model varies n, I, and T simultaneously, similar optical depths through the slab can be obtained 
for varying combinations of the parameters. However, we are more concerned with their product, 
r, than with the specific values of the parameters themselves. This interrelation between model 
parameters is discussed more fully in Appendix B (see Tables Bl and B2 and Figures Bl and B2). 
A detailed discussion of our fitting procedure and the assumptions included in our model are also 
given in the appendix. 



4.2. Main Sequence vs. Pre main Sequence Templates 



It is well known that WTTS, though not accreting, still have active chromospheres capabl e 
of producing small amounts of excess emission at UV wavelengths (e.g. iHoudebine et al.l ll996), 
Ingleby et al.1 (120111 ) recently pointed out the advantage of using WTTS templates in fits to CTTS 
spectra in order to determine more accurate accretion excesses; use of MS templates can result 
in an over-estimate of the accretion luminosity, L acc . In order to test the impact of using MS 
templates for our sample, we have modeled 8 stars using the WTTS template V830 Tau. V830 Tau 
is a K7 star so all of the modeled stars in our sample are those that required a K7 MS template 
in the original fits. Before using V830 Tau as a template it was de-reddened and convolved with 
the instrumental resolution of the McDonald and Kitt Peak intsruments to match the spectral 
resolution of our observations. A comparison of the results is shown in Table 4. 

It is obvious that the parameters generated by the two separate fits are not exactly the same. 
Overall, however, there is general agreement. The most important result pertains to the accre- 
tion rate estimates, for which both template stars produce similar values. Higher M values from 
MS templates are e xpected based on the lack of any excess emission from active chromospheres 



(jlngleby et al 



2011 



This behavior is weakly reflected in our comparison where the MS fits pro- 
duce M's that are larger on average by a factor of 1.13. In any case, the small differences found 
here do not affect the comparisons discussed in §5. We conclude that, for the purposes of this 
study, fitting with MS templates as opposed to PMS templates results in negligible changes to the 
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relevant parameters. 



4.3. Determining the Stellar and Accretion Parameters 



Once -R* is determined from the model, the luminosity of the star can be determined by 
adopting the effective temperature corresponding to the spectral type of the best-fit template star 
using the spectral type-T e fj calibrations discussed in Appendix A. The stellar luminosity and 
effective temperature can then be used to estimate the stellar mass and age by placing the star 
on an HR diagram and comparing its position with pre— m ain sequence evolutionary tracks (e.g. 



Siess et al 



2000 



Palla k Stahlei 



masses determined u sing the 



determined from the Baraffe et al 



1999; 



Siess et al 



Bara ffe et al.lll998T ). We fi nd good agreement bet ween the 



(|2000h tracks and those of IPalla k Stahleri (|l999T ). Masses 



()1998l ) tracks typically differ from the other two by 50%. This is 



not unexpected due to differences in the way each set of models treats the stellar inter ior equation 



of state and convective en ergy transport (iSiess et al.l 120001 : iHillenbrand k White! 12004 ) . The ages, 
however, determined using baraffe et all lll99sh tracks ar e systematically h igher by a factor of ~2- 



10 than those found using the lPalla k Stahleri (119991 ) and lSiess et al.1 (120001 ) models. The large age 
difference produced by using different sets of pre-main sequence eveolutionary models is important 
to keep in mind when comparing cluster ages determined by different studies. 



Hillenbrand k White) (|2004l ) find that PMS models do not agree well with dynamical mass 
estimates for M*<0.5 M and that the models underpredict actual masses by 10%-30% for 0.3 
M0<M*<1.2 Mq. Again, the uncertainty is mainly due to our relatively poor knowledge of con- 
vective energy transport and the necessary opacities at low effective temperatures. A large portion 
of our sample (80%) have estimated masses < 1.2 M , while ~33% have masses < 0.5 M Q . At 
these low masses the models become even more uncertain, due in part to the uncertain location 
of the stell a r birt hline (jBaraffe et al.ll2002l ). For the remainder of this paper, we choose to use the 
Siess et al.l (|2000l ) models in order to det ermine mass and age e stimates for our sample, though 
use of the quantities determined from the lPalla k Stahleri ([1999) models produces only small dif- 
ferences (<10%) in the final stellar and accretion parameter values (e.g. R*, M, and f acc )- More 
specifically, we use the Z=0.01 (convective overshooting not included) tracks for our mass and age 
estimates. Uames et al.1 (|2006l ) found some evidence that young clusters tend to be metal-poor, 
although the number of clusters examined in their study (3) is small. Thu s , we a dopt the slightly 
less-than solar metallicity tracks for our analysis. In any case, iMayne et ()2007l ) have shown that 
age estimates for young clusters do not diff er significan t ly bas ed on the adopted metallicity in the 
evolutionary model and a comparison of the lSiess et al.l (|2000l ) Z=0.01 and Z=0.02 (approximately 
solar metallicity) tracks shows that differences in mass estimates are negligible at comparable values 
of and T e fj. 

Figure 2 shows our sample plotted on the HR diagram. Tables 5, 6, and 7 list M*, the age, 
R*, T e jf, L*, the rotation period P ro t, face. Ay, and the derived accretion rate, M, for each star. 
Accretion rates and filling factors are not calculated for the NA stars and so are not included in 
Table 7. Uncertainties in the mass and age estimates, apart from the inherent uncertainties in PMS 
evolutionary models, are dominated by uncertainties in the underlying effective temperature of the 
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photosphere and also by uncertainties in the stellar radius, which is determined from the model. 
Typical error bars for L* and T e ff are shown on Fig. 2. The uncertainties in T e ff and L* translate 
into 30-50% uncertainties in the stellar ages and 20-25% uncertainties in M*. We also note that 
the radiative core masses used in Sections 5 and 6 are subject to similar, if not larger, uncertainties 
as M*. 

Good fits to the stellar continuum should provide reliable estimates of the underlying contri- 
bution of the photosphere. If the underlying stellar contribution is known, and T e ft is estimated 
with some confidence, the radius can be determined with typical uncertainties of ctr=±15%. Once 
R* and M* have been determined, M can be directly calculated using the following relationship 



(|Gullbring et all 11998 ): 



M =^wV-w r ) ■ (4) 

To determine M we take the disk truncation radius, Rtt, to be 5i?*. The accretion luminosity, 
L acc , is determined by adding up the slab flux F acc across the wavelength range 100 < A < 40000 
A and applying the relationship 



4"^^?^ face-Face 



(5) 



which only considers the luminosity of the accretion columns onto the stellar surface. The flux, 
F acc , is primarily determined by the observations, i.e. from the strength of the Balmer jump and 
the shape of the Paschen and Balmer continuums. 

It can be seen from Figure 2 that the SA stars, which have an average age of ~2.0 Myrs, appear 
generally younger than the WA and NA stars. The WA stars (~10.6 Myrs) are on-average similar 
in age to the NA stars (~9.2 Myrs). The NAs show a la rge spread in age and rotational period . 
Large age spreads in star fo rming regions are common (e.g. iDahm &: Simonll2005l : lHillenbrandlll997l : 
Kenvon &: Hartmannlll995l ) and the approximate age of the cluster is usually taken to be the most 
common value in the distribution. Due to the small number of stars in the NA group, we cannot 
make any definitive statements about their evolutionary status. We can, however, generalize and 
say that the older, shorter period NA stars have probably decoupled from their disks and begun to 
spin up as they contract. L4511 (~7 Mrys) and L6024 (~1 Myr), both longer period stars, have 
most likely decoupled from their disks only recently and have not had time to spin-up appreciably. 
L4098, which is both young and fast-rotating, seems to have not had time to establish disk-locking 
and thus is still contracting without being regulated by its disk. 

The clear separation in age of the SA and WA/NA stars hints at a real age diffe rence between 



the g roups. This same behavior is seen in the Taurus-Auriga star forming region: iBertout et al 



(j2007h . using revised distance estimates for 36 WTTSs and 30 CTTSs, show that the CTTSs 
are systematically younger as a group than the WT TS population a nd th us the former are most 
likely the predecessors of the latter. The result from IBertout et al.1 1)20071 ) is the first convincing 
observational evidence that CTTS and WTTS are not coevolutionary but rather WTTS are the 
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evolved products of the CTTS phase. Our HR diagram of NGC 2264 (Fig. 2) seems to su pport this 
scenario, although our sample size is approximately half that of the iBertout et al.l (|2007l ) study. 



Our sample displays a large spread in age, with some stars seemingly as old as 50 Myrs. The 
average age is ~6.4 Myrs which does not agree very well with previous, more statistically significant 
age determinations of the cluster (Mayne et al. 2007-3 Myrs; Dahm & Simon 2005-1.1 Myrs for the 
TTS population; Rebull et al. 2002-3 Myrs; Park et al. 2000-3. 2 Myrs; Sung et al. 1997-3 Myrs; 
Walker 1977-3 Myrs) . To our knowled ge, our study and that of iMavne et all (|2007l ) are the only 
ones that employ the lSiess et all (2000) models to derive age estimates for targets in NGC 2264jj 
so it is important to reiterate that different PMS evolutionary models can produce very different 
age estimates f o r the same object. As we noted above, com parisons between ages determined from 
the lSiess et all (|2000j) models and the lBaraffe et al .1 (|1998l ) PMS tracks yield values that differ by 
up to an order of magnitude; individual ages should be examined with caution. 



4.4. L H p vs. L acc 

We have measured the H/3 line luminosities for our sample. The line flux was calculated by 
subtracting the continuum from 4820-4840 A averaged with the continuum from 4880-4900 A and 
summing up the residuals from 4840-4880 A. A distance of 760 pc is assumed when computing the 
line luminosities (see Appendix B for a discussion of distance estimates) . The Ay values from Tables 
5-7 were used to correct each star for extinction. We used our SA and WA stars to determined the 
relationship between the H/3 line luminosity and the accretion luminosity. Figure 3 shows a strong 
relationship between the accretion luminosity, L acc , and the H/3 line luminosity, Ljjp- The best-fit 
relationship to the data is 



log 



(2.80 ±.21) + (1.12 ±.06) log 



L 



Hp 



L, 







(6) 



Equation (6) is almost identical to the relationship calculated bv lFang et all (|2009l ) for TTSs 
in the Taurus- Auriga star forming region. Eq. (6), however, was calculated by excluding any stars 
with logio (Lh@ / Lq) < -4.50 in order to account for chromospheric line emission not produced in 
the accretion flow. This number was chosen by averaging the H/3 line luminosities of 5 WTTSs in 
Taurus (Table 8), each of which show no evidence of a circumstellar disk at near-IR or millimeter 
wavelengths and therefore are most likely experiencing negligible accretion. Below this threshold 
the H/3 luminosity is likely dominated by chromospheric emission. Therefore, it is not possible 
to measure the H/3 luminosity produced by accretion for TTSs that lie below this value and are 
within the mass and radius range studied here, i.e. H/3 cannot be used as an accretion proxy for 
these stars. The NGC 2264 objects with line luminosities below this value are not included in the 
determination of eq. (6). Also plotted in Fig. 3 are the NA stars (red diamonds) with measureable 
H/3 emission. They are plotted using the accretion luminosities that result from fitting the spectrum 



Rebull et all (12 002) 



comp are the ages obtained for NGC 2264 using the Siess et al. 1 2000l l models with those of 
D'Antona fc Mazzitellil (1994) but they choose to adopt values from the latter models. 
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with a template and a slab. All of these objects are clustered near the threshold and lie below the 
line of best fit (eq. 6), thus indicating that H/3 line strength should not be used as the sole proxy 
for the accretion luminosity. The NA stars are not included when determining eq. (6). 



The Fang et al. (j2009l ) relationship includes ~7 objects below our estimated threshold. We 
do not know the identities of these object jf] but they are most likely late M-type stars and brown 
dwarfs from lHerczeg fc Hillenbrand! (|2008l ) in which case our threshold estimate may not apply if 
the chromospheric contribution to H/3 decreases at very low masses. The application of eq. (6), 
however, to a large sample of TTSs in the same mass range as our sample and with Lnp less 
than the threshold value is cautioned against due to the intrinsic H/3 luminosity produced by active 
chromospheres. Doing so can lead to unrealistic accretion luminosity estimates and in turn to 
incorrect mass accretion rates. 



4.5. Comparison with Rebull et al. (2002) Values 



Rebull et al.1 (120021 ) (hereafter R02) computed the age, M*, R*, Ay, 



and M for 160 of the 36 

stars in our sample. A comparison of the values is given in Table 9. Visual comparisons of Ay and 
M are plotted in Figures 4 and 5, respectively. Differences in the derived masses range up to 75% 
but, on average, they agree to within 30%. The most striking difference between the parameters 
is the age determinations. The large di s parity can be primarily attributed to the us e of different 
PMS evolo tionary models. iRebull et all (120021 ) use the iD'Antona h Mazzitellil (119941') tracks while 



we us e the lSiess et al.1 (|2000h tracks. As |R02j point out (IROl Fig. 13) . the 



D'Antona & Mazzitelli 



( 1994 ) models consistently produce ages younger than the Siess et al. (j200ol ) models by .5-dex. 



The Ay values derived by |R02| for several stars using (R-I) color excesses are significantly 
different than our Ay values, which are the best-fit value s fro m our models (Figure 4 ). W e have 



R02 value, 



attempted to model our stars using the Ay values from |R02| . If Ay is fixed at the 
however, the fit is either not of comparable quality to the fit determined by letting Ay vary or 
our routine fails to find an acceptable match to the data. For a few of th e sta rs, the difference 
in Ay can be attributed to differences in the spectral type determined by IR02I and the spectral 
type of the template star used to model the spectrum. In order to derive accretion rate estimates, 
R02 use E/-band mag nitudes to estimate the accretion luminosity base d on the relationship from 
GuUbring et al.l ()1998l ). They are well correlated (Figure 5), with |R02j producing accretion rates 
higher on average by a factor of ~2.5, though the factor is larger for stars with smaller accretion 
rates and smaller for stars with higher M. 



There is one discrepency, however, that deserves closer attention. IRebull et al.l (|2002l ) calculate 
a large accretion rate (M=1.3 x 10 Mq yr _1 ) for the star L6172, while we determine, based on 
L6172's optical spectrum, that it is not significantly accreting. Using only a scaled and reddened 
template star, we find an excellent fit to the L6172 data. If we allow some accretion luminosity 



^Fang et alj J2009I ) take their H/3 luminosities from the studies of Herczeg fc Hillenbrand ( 2008I ) and 
GuUbring et"aD jl998l ). 



Values of Av and M are not calculated by R02 for a few individual stars in this subsample. 
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to contribute to the spectrum, fits of comparable quality can be obtained which produce very low 
accretion rates of < 8 x 10 -10 M yr -1 . A glance at Fig. lc confirms the lack of any obvious excess 
emission redward of 3400 A. Fig. lc also shows that L6172 has reasonably strong H/3 line emission, 
but this could be du e alm ost entirely to chromospheric activity. Using eq. (6) and the values of 
M*, i?*, and M from |R02| . we have plotted L6172 in Fig. 3 (open blue circle). The H/3 luminosity 
calculated using eq. (6) is logio(Lff/3/L )=-2.59, while we find logmiLnp/ Lq)=-3.76, more than 
an order of magnitude less than the value calculated using the IR02I values. Although TTSs are 
known to show significant variability on a variety of timescales ( Bertout 19891 ). routine variability 
studies typically do not fi nd variations in the accretion rate of 3 orders of magnitude. For example, 
Alencar &; Batalhal ((2002) find accretion rate variations for TW Hya of ~1 orde r of magnitude over 
a year timescale. In a UV variability study of BP Tau, lArdila Basril (|2000l ) measure accretion 
rates that vary by a factor of 7 over ~1 year timescales. The FU Orionis phenomena can produce 
increases in the mass accre tion rate by 4 orders of magnitude or m ore but the decay timescale is 
often longer than a dec ade (lHerbiall977l ; lHartmann fc Kenvonlll996l ). Our observations were taken 
~8 years after those of |RQ2j so the possibility of L6172 having undergone an intense period of mass 
accretion during the latter set of observations cannot be entirely ruled out. 

Although it is not clear to us what the source of th is discrepency is, based on the relative 
agreement of our other accretion rates with those of IRQ 21 . we have confidence that L6172 is not a 
heavy accretor in our data for 2 reasons: 1. It does not appear to have any excess emission in the 
optical, and 2. P ro t for L6172 is 1.76 days which implies that it is most likely rotating freely and is 
not locked to its disk. Although L6172 appears to be very young {tl§\72 < 2 Myr in both studies), 
which makes it more unlikely that it has ceased interacting with its disk, it is certainly possible 
that it has yet to establish disk-locking, or is in the proc ess of doing so. This is consistent with 
the "moderate" angular momentum loss scenario used by lL05l to explain the presence of young, 
fast-rotators with large IR excesses in Orion and NGC 2264. We believe that L6172 falls into this 
category. Aside from this individual case, we find good agreement within the (large) uncertainties 
between our parameters and those of IR02 . 



4.6. A Note on Uncertainties in M 

We have already discussed the uncertainties in our model fits (see Appendix B for further 
discussion). Calculated values of M rely on the stellar and accretion parameters derived from 
the models. Assumptions concerning the geometry of the system and the adopted value for Rx r 
introduce a factor of ~2 uncertainty in individual accretion rates. For example, we assume that the 
observer is looking straight "down" on the accretion column instead of modeling multiple accretion 
columns at different angles and projections onto the surface. We also assume that half of the 
accretion luminosity is emitted into the solid angle subtended by the star, while half is emitted 
outward toward the observer. These assumptions affect the final values of L acc and f acc which 
translate into uncertainties in M. A more obvious uncertainty arises in the choice of Rt t which 
can result in differences in M of a factor of ~1.6 depending on the adopted value in the range 3-6 
i?*. In addition to these uncertainties, we can estimate uncertainties in M, within the context of 
our assumptions concerning the accretion flow and system geometry, once we estimate uncertainties 
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in the derived parameters. 

Our x 2 fitting routine estimates uncertainties in the model parameters based on the quality 
of the fit. We can easily translate these uncertainties into uncertainties for R* and M*. If we 
assume an uncertainty of 150 K in T e ff, uncertainties in M* are generally near 20-25%. Typical 
uncertainties for R* are 15%. For accreting stars, L acc is relatively well constrained by shape of the 
spectrum and its measured strength, although i?* and f acc are used to calculate the final value. To 
compute the uncertainty in L acc , we assume that the uncertainty in F acc is equal to the errors in 
the observed spectrum. 

It is important to highlight the role that Ay, the visual extinction, plays in determining the 
accretion rate. Higher values of Ay will cause the bluer portion of the spectrum to be enhanced 
relative to the red wavelengths thus essentially increasing the slab luminosity required to generate a 
good fit to the spectrum. A higher accretion luminosity will result in a higher value of M assuming 
that R* and M* remain the same. As we demonstrate in Appendix B, we are relatively confident 
in the best-fit parameters for each spectrum, i.e. the spread in the final values of each parameter 
is small when the input parameters are varied. This is true for Ay, as can be seen in Table A2. 
As m entioned in §4.5, however, our model Ay values can differ significantly from those derived by 



R02l (Fi gure 4). The disagreement between our model Ay values and those determined directly 



from the R02 photometry are not surprising considering th e different wa v eleng th regions utilized in 
the calculations. As we discussed in §4.5, fixing Ay at the iRebull et al.1 (120021 ) values when fitting 
the spectra generally does not produce models of the quality obtained by letting Ay vary. We 
highlight these points simply to illustrate the caution necessary in considering any individual value 
as a precisely determined quantity. 

Combining all of the uncertainty estimates yields uncertainties in M of 35-50%. Thus, uncer- 
tainties in M produced by uncertainties in the stellar and accretion parameters are small relative 
to the uncertainties in the treatment of the accretion flow itself. For this reason, individual values 
of M should be view with caution. However, most of the systematic uncertainties, as mentioned 
previously, should not affect the final correlations presented in Section 5: that is, such uncertainties 
may shift the stars as a group, but it is unlikely that the errors that may be present will introduce 
a correlation when none is actually present. 



5. ANALYSIS AND RESULTS 

In order to test the magnetospheric accretion models outlined in §2, we use the stellar and 
accretion parameters from Tables 5 and 6 to plot the left-hand side vs. the right-hand side of the 
relationships predicted by eqs. (1) and (3). We then test whether or not the data are correlated. 
The NA stars are not included in our analysis because they are not accreting measureable amounts 
of gas from their disks and thus are most likely no longer strongly coupled to their disks. The same 
analysis is performed for both the assumption of a constant magnetic field from star to star and 
our assignment of magnetic field strengths to each individual star based on the size of the star's 
radiative core. 

We have discussed the uncertainties in our derived parameters at some length. Most of these 
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uncertainties are systematic and dominate the observational uncertainties to a large degree. For this 
reason, we expect the measured errors to produce scatter about the predicted relationships while 
the systematic uncertainties will tend to shift the relationships but not create correlations where 
none actually exist. Scatter about the predicted relationship is probably also due to fluctuations in 
the equilibrium configuration of the CTTS system . In reality, CTTS system s most likely experience 
variable accretion rates (e.g. lArdila Basrill2000l : Ivan Boekel et al.ll2010l ) leading to variations in 
M, P ro t, and RTr, the disk truncation radius. Long- term variations in the stre ngth of the dipolar 
component of the magnetic field (e.g . V2129 Qph: seelDonati et al.ll2007l . l2011al ) may also affect the 



equilibrium star-disk configuration (jPonat i et al.ll2011al lbl) which may in turn affect the correlation 



measurements. Thus the correlations predicted by eqs. (1) and (3) are likely only valid for the 
average of these fluctuations around the equilibrium. This effect can be seen in Fig. 8c where a 
strong correlation is present and also exhibits the expected scatter surrounding the line of best-fit. 

In order to quantify the strength of the relationships given by eqs. (1) and (3) as applied 
to our sample, we compute the S pearman's rho rank correlation coefficient, r, and the correlation 
significance P (jPress et al.1 119871 ) for three separate group of stars: the strong accretors, weak 
accretors, and the combination of strong and weak accretors. These results are shown in Figures 
6 (eq. 1) and 8 (eq. 3) for th e case of a constant magnetic field for the sa mple . We also tested 
eq. (1) using stars from lR02j for which rotation periods are available from |L04| . This sample is 
shown in Fig. 7. In order to prevent extreme values from overly influencing the correlations, we 
consider the logarithm of the quantities in the proportionalities. The correlation significance, P, 
is a measure of how likely it is that the calculated correlation is not actually present in the data. 
Thus, a lower value for P implies a lower liklihood that the correlation is not real. Values of P<.01 
(99% confidence level) are generally considered evidence of a real correlation. Values of r close 
to -1.0 and 1.0 imply strong correlations. A summary of the correlations for each equation, their 
associated values of P, and the slope of the best-fit line are given in Table 10 for the case of a 
constant magnetic field strength from star to star. Best-fit lines are not calculated for the poor 
correlations produced using eq. (1). 

It is obvious from Table 10 and Figs. 6a, b, c , and 7 th at only weak cor relations, if any, 
exist for eq. (1), which is derived from the iKonigll (119911 ) and I Shu et al.l (119941 ) magnetospheric 
accretion models by assuming a purely dipolar magnetic field and a stellar rotation rate equal to 
the Keplerian rate in the disk at Rt t - It is equally clear from Figs. 8a, b, and c that eq. (3), 
a relationship predicted by |JG02r s modified version of the IOS95I theory to include a non-dipolar 
surface magnetic field, yields excellent correlations with low values of P for all three samples. Again, 
we have assumed a constant magnetic field for our sample when plotting the proportionalities in 
Figs. 6, 7, and 8. 

We applied the same analysis to our sample after assigning magnetic field strengths to each 
sta r based on th e extent of the star's radiative core, following the suggested evolution o f the dipole 



bv iDonati et al.1 (12011al ) . The internal structure of each star was estimated using the ISiess et al 



20001) models. Usin g t he cu r rent o bservational results of dipolar field strengths on TTSs as a guide 



Donati et al.l 1201 lal Jbl. l2010i . 120081 ). we chose values of 0.1, 0.5, and 1.5 kG to be assigned to stars 
with M core /M*>0A, 0.0 < M core /M*< 0.4, and M core / M*=0.0, respectively, where a star with 
Mcore/M* =0.0 is fully convective. In other words, instead of assuming a constant magnetic field 
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from star-to-star we let it vary according to its internal structure to see the overall effect on the 
correlations. We note that the values given above for the surface dipolar field are values at the 
magnetic pole, while eq. (1) uses the equatorial field strength. For the dipolar field component, 
the equatorial value is one-half of the polar value. Although eq. (3) involves the surface magnetic 
field value participating in the accretion flow, independent of its specific geometry, we have used 
the same field values given above as for the purely dipolar case of eq. (1). After letting the 
magnetic field values vary, the correlations actually worsen for both equations: for eq. (1), we 
now have r=0.13 and P=0.52; for eq. (3), we now have r=0.75 and P=2xl0~ 5 . The variation of 
the magnetic field essentially spreads out the data along the abscissa and makes the best-fit slope 
shallower (slope=0.50) so it is not surprising that the correlations become weaker. Thus, while the 
correlation strengths change when the dipole component of the magne tic field is allowed to vary 
from star-to-star, our results still show stronger support for the IOS95I model than for the purely 
dipolar magnetospheric accretion models from eq. (1). For the remainder of this paper we will refer 
to the constant magnetic field analysis for both eqs. (1) and (3) when discussing the implications 
of our results. 

The slope of the predicted relationship should be 1.0 using eqs. (1) and (3). For the initial 
linear fit to the data, we assume that each point has an equal uncertainty. The standard deviation 
of the residuals between the best-fit line and the data are then assumed to be the la uncertainties 
for each point. These la values are then used to compute the uncertainty in the slope itself. As 
can be seen in Figure 8c and Table 10, the slope (0.63±0.03) of the best-fit line to the data (solid 
line) does not match the predicted line (dashed) slope of 1, even to within 3a. IJG02I also find a 
shallower slope (0.73=L0.14) for their sample of stars in Taurus. 

The cause of the slope differing from the predicted value in Taurus and NGC 2264 is currently 
unkown; however, one possibil ity may have to do w ith the evolution of the s tellar magne t ic fiel d 
geometry. As described in §2, iDonati et al.l (|2011bl ). building off the work of iMorin et al.l (|2008l ). 
suggest that very young, fully convective PMS stars have field geometries with relatively strong 
dipole components. As they age and develop even a small radiative core, Donati et al. suggest 
the surface field becomes much more complex, with higher order terms becoming stronger than the 
dipole. As a result, younger stars will have magnetic interactions with th eir disks at several stellar 
radii that are essentially the same as for a pure dipole as assumed by lJG02l . However, for older stars, 
the higher order magnetic components may start to change the nature of the field interaction with 
the disk and thereby change the specific predictions of how stellar and accretion parameters should 
scale with one another. This could result in a changing of the slope of the predicted relationships. 
If true, the youngest clusters should show correlations with slopes close to 1.0 and t he slop e should 
evolve away from this value for older clusters. This is essentially what is seen with |JG02| finding a 
slope closer to 1.0 for their Taurus sample than we find for our slightly older NGC 2264 sample. 



At first glance, our data appears to strongly support lJG02l 's modified version of the IOS95 
theory . Before we are able to interpret the strong correlations as conclusive support for the modified 
OS95I model, however, we must explore any inter dependencies of the variables that may arise for 
other reasons. In this way, we can ensure that the observed correlations are in fact evidence for the 
model predictions and are not driven by expected relationships between parameters based on how 
they are determined. 
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5.1. Testing the Correlations 

For a uniform flow, the mass accretion rate is 

M = A7rR 2 J acc pv (7) 

where v is the velocity of the flow and p is the density, both measured at the point of impact on 
the stellar surface. It is clear that M should be correlated with R^f acc - Thus it is important to test 
whether or not M and i?^, M and faccR*, and M and face, show any signs of being correlated in our 
sample. Table 11 shows the results of these comparisons. Clearly there is a significant correlation 
between M and R^ (Figure 9) and a lack of correlation between M and f acc . In fact, one might 
worry that t he corr elation shown in Fig. 8c for the full sample of SA and WA stars plotted using 



the modified IOS95I prediction is dominated by the correlation between M and the product of f acc 
and R*. However, by comparing the correlations from Table 10 and Table 11 we can see that the 
inclusion of M* and P ro t does enhance the quality of the relationship when looking at eq. (3). 
The expectation of a correlation between M and R^ is obvious from eq. (7), but the lack of a 
correlation between M and f acc is not quite so clear. Furthermore, if the good correlation seen in 
Figure 8 is driven primarily by the correlation between M and R%, we would also expect this strong 
correlation to show up in the plots for eq. (1) (Figure 6). As we have seen, however, there are no 
correlations present in the eq. (1) plots leading us to the conclusion that the M—R% correlation is 
not the primary driver for the correlations for eq. (3). 



A further test of eq. (3) is given by IJG02I as 



R* FaccProt ( 8 ) 

which can be obtained by substituting eq. (4) into eq. (3), solving eq. (5) f or f nr( - and plugging 



the result into eq. (3). Equation 8 provides the best test of the modified IOS95I theory in the 
sense that the parameters involved are determined as independently as possible from one another. 
Although both L acc and F acc are ultimately determined by model fits to the data, F acc is the 
quantity directly measured from the observations. In addition, L acc and F acc are not well correlated 
(r=0.21, P=0.31) so any real correlation shown by eq. (8) is a result of the full comparison of all 
the combined variables. Figure 10 shows a scatter plot of eq. (8). While there is some evidence 
of a correlation (r=.46), the probability of no correlation is high enough (P=0.02) to prevent any 
definitive statements about the support of eq. (3), rewritten here in the form of eq. (8). However, 
when L7974 is removed (the lower right-most blue circle) from the analysis, the correlation improves 
to r=0.65 with P=.006 suggesting a stronger relationship between the variables. 

A similar test can be done for eq. (1). By replacing M in eq. (1) with eq. (4) and dropping 
the constants, eq. (1) can be rewritten as 

RT oc MlHl&PlS. (9) 



Figure 11 shows the results of plotting the SA and WA stars using eq. (9). There is no correlation 
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present at all, confirming the lack of support for the dipole models of eq. (1). In fact, the correlation 
is worse than the correlation from eq. (1) (Figure 6) indicating that the correlation between M 
and R% is most likely bolstering the weak correlation found for eq. (1). 

We have also tested the impact on the correlations of calculating individual values of Rt t for 
each star. Each truncation radius was computed by assuming the star is locked to its disk. In other 
words, Rtt is simply the radius at which the disk rotates with the same Keplerian rotation rate as 
the star. The changes in M using the individual values of Rt t range from as little as 1% to as much 
as 50%. However, when the analysis is repeated using the new values of M, the correlations do 
not show any significant change for either eq. (1) or eq. (3) (eq. l-r=0.33, P=0.10; eq. 3-r=0.90, 
P=7xl0~ 10 ). This is not surprising since changing the truncation radius by even a factor of 10 will 
most often result in less than a factor of 2 change in the accretion rate. 

We believe that the modera te correlation shown by eq. (8) provides additional support for 
the modified version of the IOS95I theory of magnetospheric accretio n. The very strong correlations 
shown in Fig. 8, however, provide the strongest support for the IOS95I model even though we 
find that the eq. (3) correlations could possibly be driven by a correlation between M and R^. 
If the correlation did not improve with the inclusion of M*, P ro t, and fg nr . we would hesitate to 
draw any firm conclusions concerning the validity of the modified IOS95I theory. However, since 
the correlations do improve and the result is a very tightly constrained relationship between the 
two sides of eq. (3), it appears the analys is of |JG02I to include non-dipole surface magnetic field 
geometries within the context of the IOS95I model is supported by the data in NGC 2264. 



6. DISCUSSION 



Using stellar and accretion parameters for TTSs in the T aurus-Auriga star forming region 



JG02I also find support for the modified version of the OS95 theory through the same type of 
analysis we present here. The correlation presen ted by I JG02I for eq. (8) (their eq. 12) is stronger 



than the one we find for our sample. The UG02I correlation appears to not be strongly influenced 
by a correlation between R% and M, while the effect of such a relationship may have a larger 
influence on ou r results for eq. (3). Nonetheless, bot h studies show much stronger support for th e 
mo dified OS95I theor y than they do for the models of Konigl (1991), Cameron Sz Campbell ( 1993 ). 
and Shu et al. all of which assume a purely dipolar stellar magnetic field. As IJG02I point 

out, this does not rule out the domination of the dipolar component at radii comparable to Rt t due 
to the faster decrease in field strength of higher-order components with increasing distance from 
the stellar surface. In fact, S. G. Gregory et al. (2012, in press) show that the dipole component 
dominates at the truncation radius even if the field strength of the octopolar component is ten 
times that of the dipolar component at the magnetic pole. As mentione d briefly in §2, there is 



mounting evidence that magnetic fields on TTSs are not purely dipolar (jValenti &: Johns-Krull 



2004 



Johns-Krull 


2007; 


Donati et al. 


2011a 


b 



b c])- While higher order magnet i c field components 



2012J), estimates of the disk truncation radius can be made using the dipolar component only and 



ignoring the negligible contributions from the higher order components. 
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Based on their rec ent measur e ments of magnetic fields for a small sample of TTSs (e.g. 



Donati et al. 


2011c 


b 


a, 


2010. 


2008, 


2007) 



201(3 . 120081 . 120071 ). Donati et al. have speculated that the dipole com- 



ponent of the magnetic field should become weaker relative to the higher order field components as 
the star ceases to be fully convective and develops a radiative core. This evolution of the magnetic 
field geomet ry is given as eviden ce that the magnetic fields of TTSs are dynamo-driven rather than 
primordial (IDonati et al.ll201ld ). If the strength of a star's dipolar magnetic field decreases, it is 
expected that the accretion disk would be able to move in to a smaller radius enabling the star 
to spin-up and decrease its rotation period. Thus, according to the suggestion of Donati et al., 
one might expect the rotation periods of partially convective stars to be shorter as a group than 
the fully convective stars. Eq. (1), however, shows that the relationship between B^ip and P ro t 
also depends on the star's mass, radius, and accretion parameters, i.e. suggesting that P ro t should 
decrease as B^ip decreases is an oversimplification of the problem. Therefore, in order to test this 
hypothesis, we have plotted (F ig. 12, upper panel ) the s trength of the equatorial dipolar magnetic 
field component as derived by iJohns-Krull et al.l (|1999l ) for the IShu et al.l (|1994l ) magnetospheric 
accretion model (effectively solving eq. (1) for B^ p but with the appropriate constants to turn the 
proportionality into an equality) against the ratio of the mass of the radiative core to the total 
stellar mass. This is the magnetic field required by disk-locking theory to maintain the observed 
rotation rate given the measured m a ss, ra dius, and accretion rate onto the star. The core mass 
was calculated using the ISiess et al.l (|2000i ) evolutionary models for the SA and WA stars, where 
M core /M* =0.0 is a fully convective star. The NA stars are not included due to the lack of derived 
accretion rates required to estimate B dip- 
It can be seen from the upper panel of Fig. 12 that the derived dipolar component of the 
magnetic field at the stellar equator does not correlate with the size of the star's radiative core. 
We have also examined the predicted field strengths of the fully convective and partially convective 
samples using a Kolmogorov-Smirnov test. The significance of the K-S test is Pr:5=0.24, large 
enough to prevent the rejection of the null hypothesis that the two samples are drawn from the 
same distribution. We also performed a K-S test on the rotation periods of the two samples, which 
are shown in the lower panel of Fig. 12. In this case, Pa:5=0.38, again providing little evidence that 
the samples are drawn from separate distributions. The upper panel of Fig. 12 was made assuming 
two things: 1. the star-disk interaction is dominated by the dipolar magnetic field component, 
and 2. the star is locked into co-rotation with its disk, i.e. disk-locking. Due to the lack of 
correlation seen in the upper panel of Fig. 12, one or both of these assumptions must be incorrect 
or the evolution of the magnetic field as suggested by Donati et al. is incomplete. However, our 
data show good support for disk-locking and the dominat ion of the star-disk interaction by the 



dipolar magenetic field comp onent is well supported (e.g. lAdams k, Gregory 



2012: 



Donati et al 



2011a: 



Mohantv fc Shul l2008). More measurements of the dipolar field component are needed for a 



larger sample of TTSs in order to reinforce or reject the recent findings of Donati et al. 



While our data (and that of IJG02l ) support a somewhat modified version of the IOS95I model 
and therefore indirectly provide suppo rt for the magnetocent rifugally driven wind (X-wind) which 
is a component of that model (see also lMohanty &: Shull2008l ). other competing theories have been 
developed to explain angular momentum loss in low mass stars and thus the observed slow rotation 
in most CTTSs. These models invoke an accretion driven stellar wind as the main momentum 
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regulation device (e.g. iMatt Pudrita 120051 . 12008 ) and argue, based on numerical calculations of 
magnetic torque interactions between the disk and star, that the torques associated with disk- 
locking theories are not sufficient to produce the observed slow rotation rates. In fact, when only 
considering t he magnetic coupl ing of the star and disk, these theories predict a positive torque on 
the star (e.g. IMatt et al.ll2010l ). causing it to spin up. Using the same analysis discussed in § 5 
we have tested the equilibrium equations of lMatt &; Pudrita (|2008l ) (Matt, private communication) 
involving M*, M, R*, and P ro t- For the case R rn ^R,T r, where R r „ is the corotation radius, we find 



weak evidence of a correlation (r=0.42, P=0.04). The lMatt fc Pudrita (|2008l ) relationships assume 
a purely dipolar magnetic field. Thus, based on the lack of support for eq. ( 1), it is not surprising 
that only weak correlations are present in our data for the lMatt &: Pudrita (|2008l ) relationships. 



We beli eve ou r analysi s provides suppo r t for the IOS95I model of magnetospheric accretion, 
extended by lJG02l . and later iMohanty Shul (120081 ). to include non-dipole field geometries. These 
results also support the disk-locking scenario which see ms to be well supported observation ally 
based on period distributions of young stellar clusters (e.g. 
thou gh conflictin g observations do e x ist (e.g. 



Stassun et al.1 



erbst et al. 2002 



1999; 



Lamm et al. 2005), al 



Rebull et al.1 



. 2002; 



Makidon et al 



20041 ) . However, iRebull et al.1 feood ) , Icieza fc Baliben faooj ) , and bahm. Slesnick. fc Whitel feoij ) 



find the strongest evidence to date that circumstellar disks are directly involved in regulating the 
angular momentum of their central stars, effectively ruling out a mechanism that does not involve 
significant angular momentum exchange between the star and its disk. Accretion-driven stellar 
winds could certainly provide a m eans o f angular momentum removal for CTTSs but the signifi- 
cant correlations found here and by lJG02l for eq. (3) provide observational support for the existence 
of disk-locked systems, a constraint that is not necessary in the stellar wind models. 



7. SUMMARY AND CONCLUSIONS 



We have derived stellar and accretion parameters for 36 TTSs in NGC 2264 using spectropho- 
tometry measurements taken in 2004 and 2005. Our estimate for t he age of the sam ple (~6.4 
Myrs), calculated using the pre-main sequence evolutionary tracks of ISiess et al.l (|2000l ). is older 
than several more statistically significant age determinations which aver age to ^3 Myrs. This is 



due in large part to: 1. our small sample size, and 2. our use of the ISiess et al 



model s whi ch tend to produce older ages c o mpared with the commonly used tracks of 



fll99J ) and iD'Antona fc Mazzitellil (11994! ) (jEebull et al 



2002 



accretion rates for our stars are similar to estimates by 



R02 



Dahm Simonl 120051 ) . The mass 



(200 0|) PMS 



Baraffe et al 



though the agreement is worse at 



lower M. Using t he derived pa ramet ers we test analy tic predictions from the purely dipolar disk 



lo cking m odels of iKonigll (ll99ll ) and 



of JG02's modified version of the 



Shu et al 



(1994) (eq. 1) and the non-dipolar field prediction 



Ostriker k, Shul (|1995l ) model (eq. 3). We find good support for 



the modified lOstriker Shul (|l995l ) model of magnetospheric accretion and disk-locking, although 
the correlation is influenced to some degree by a strong relationship between M and R^. A lack of 
support for the dipolar theories, however, highlights the need for an extra constraint in the theory 
or the abandonment of disk-locking. This constraint is provided by the inclusion of f acc in eq. (3). 



Although the support we find for the lOstriker Shul ()1995l ) theory is not without uncertainty. 
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our results confirm the findings of |JG02l and provide more evidence for disk-locking than against it. 



In addition, our results find no support for theories that assume a dipolar magnetic field geometry 



Daou. Johns-Krull. & Valenti 


2006; 


Mohantv k Shu 


2008; 


Donati et al. 


2008, 


2011b) 



This scenario does not exclude the possibility of accretion powered s t ellar winds as t he agent which 



removes angular momentum from TTSs (e.g. iMatt Pudritz 



2005 



Cranmerll2009l ). In fact, both 



stellar winds and disk winds launched from nea r the truncatio n poin t probably play a role in the 
removal of angular momentum, as suggested by lEdwards et al.1 (|2006l ) . 



The primary assumption made in our initial analysis is that the magnetic field strength does 
not vary significantly from star to star. We have attempted to roughly account for differing dipolar 
magnetic field strengths from star to star by assigning values based on the size of the star's radiative 
core as suggested by the recent results of Donati et al.; however, doing so actually weakened the 
correlations present in the data. In order to make this analysis more robust, better statistics on 
the strength of dipolar magnetic field components in TTSs, especially those at the ~3 Myr age, 
are needed. Whether the constant field strength assumption is justified or not for stars at similar 
evolutionary stages may help explain the slope evolution observed in our data. In addition, we do 
not find any relationship between the theoretical dipolar magnetic field component and radiative 
core mass in our sample. Future studies of emission and absorption lines from TTSs environments 
will help place constraints on the launching region of the outflows, thus helping to confirm or reject 
the hypothesis of disk-locking, and subsequent removal of angular momentum by a disk wind from 
the truncation radius, supported in this work. Until more observ ational evidence can b e gathered, 
it appears that the non-dipolar magnetspheric accretion model of lOstriker fc Shul (119951 ) remains a 
strong candidate for explaining the observed relationship between stellar and accretion parameters 
of CTTSs at the ~3 Myr evolutionary stage and, by extension, that the disk-locking scenario is 
taking place in young stars. 



A. Spectra and Model Fits 

This appendix contains the flux-calibrated spectra and model fits for our entire sample (Fig. 
Al). The SA stars are displayed first, followed by the WA and NA stars. Note that model fits 
are not given for the NA stars. The black solid line is the observed spectrum; the model fit to the 
spectrum is overplotted in red; the template star is plotted as a green solid line; the slab flux is 
plotted blue dashed line. 



B. Fitting the Data 



The best-fit models to our spectra are generated using a non-lin ear least squares fitting routine 
based on the Marquardt method (see iBevington Robinson! Il992l ). We generate our fits using a 
total of 6 free parameters: n, T, I, a, f3, and Ay where a is a flux scale factor relating the slab 
emission to the template star emission and f3 is a scale factor relating the combined slab and 
template model to the observed TTS in NGC2264. These scale factors account for the stellar 
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radius and distance (for the TTS and template) and for the filling factor of the accretion zones. 
The scale factors have the following forms: 



a 



fa 



1 



face) 
2 



R T 

dr 



'-<£)'(£)' 



(Bl) 
(B2) 



where Rt is the radius of the template, R* is the radius of the star, cLt is the distance to the template, 
and is the distance to the star. The values Rt, dr, and are known or assumed quantities. Thus, 
the final values of a and f3 yield the surface filling factor f acc and stellar radius R*, respectively. 
The template radii were taken from the Catalogue of Apparent Diameters and Absolute Radii of 
Stars (CADARS) (jFracassini et al.ll200ll ). We also calculated the template radii using published 



B-V color s, parallaxes, bolomet ric correct ions based on spectral type and effective temperatures 
taken from lcohen fc Kuh] (jl979h (K0-K61. liesselll (|l99lh (K7-M1), and lwilking. Green, k Mever 
(|1999l ) (M2-M7). All of the templates are MS stars within 30 pc. Typical uncertainties in T e ff are 
±150 K; for the bolometric corrections used, ABC=± m 05. For typical photometric uncertainties of 
5% and errors of 5% in the parallax measurements, the uncertainties in template radii calculations 
are ~15%. The uncertainty in the distance to GJ596A is much higher, at ~25%, which translates 
into an uncertainty of ~50% in the stellar radi us. By com paring the observed B-V colors of the 
templates with calibrated B-V MS colors from lCoxl ( 2000 ). we find that reddening is negligible to 
all of the templates. In most cases our calculated radii agree with the CADARS radii to within 
5%. In the cases of HD45088, GL394, and GJ596A the discrepency is closer to 15%-30%. For 
consistency, we chose to adopt the CADARS values for these templates. Table 3 lists the template 
stars and their parameters. 

The dist ance to NGC 2264 has not been p recisely determined. C urrent published values range 
from 750 pc (jMayne fc Naylorlhoosl ) to 950 pc (jFlaccomio et al.lll999h w ith a median value of ~800 
pc. The most recent determination is 910 ± 110 pc by iBaxter et al.l ((2009). Here we choose to 
adopt a distance of dNGC22&i = 760 pc. We have varied the distance used in our models from 
700 pc to 900 pc and found deviations of < 10% in the values of the stellar accretion rates, filling 
factors, and radii that result. In addition, because the distance to each individual star does not 
vary, the chosen distance to the cluster has no effect on the final parameter correlations and thus 
does not affect the main results of this study. 



The overall reddening to NGC 2264 has been estimated by multiple studies (e.g. IWalker 



1956 



Sung et al.l ll997: 



was done by 



Rebu 



Walker 



1 et al.ll2002l ). The first photometric determination of the cluster reddening 



19561 ) using main sequence O and B stars and was found to be E(U 



y)~0. m 08. lYound (jl978l ). using the standard extinction law R=3.1, determined that the cluster 



was differentially reddened and that dust clouds within the cluster varied with position, causing the 
reddening to change moderately from region to region. This was interpreted as the result of intense 
radiation from higher mass stars cle aring ou t the s pace around them and creating regions of higher 
and lower dust extinction. Indeed, lYoung) (119781 ) identifies the excessively r eddened star s in h is 
sample as those that lie behind known dusty regions of the cluster. However, iRebull et al.l (|2002l ). 
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using a much larger sample of stars, find no evidence of spatially dependent reddening within the 
cluster. This may be due in part to their use of R — I excesses which are less sensitive to foreground 
absorbers. 

Reddening values to individual stars can be affected by the local environment of the object, 
causing specific values of Ay to differ significantly from one star to the next. This individual 
reddening is what we attempt to estimate using our models. ISung et al.l (119971') confirm the dif - 



ferential nature of the reddening across the cluster found by lYound (| 19781 ) . ISung et al.l (|1997l ) 
also determine E(i? — V") to be_0. m 071±.033, in agreement with earlier estimates. In the more 



recent study of iRebull et al.l (J2002J), the cluster reddening is found to have a slightly higher value 
of E(i? — 1/)=0.146±0.03, corresponding to ^4^=0.41. These estimates, however, are "most likely" 
values and cannot be used in the context of individual stars. Here, we solve for individual model 
estimates of Ay. 

Strong emission lines from elements other than hydrogen can affect the fits to the data. For 
this reason we ignore ~10 A of spectrum on either side of any strong emission lines. The most 
common examples are the Call H & K lines at 3933 and 3969 A. The same procedure is applied to 
the Balmer lines discussed in §4.1. Ignoring these lines provides better fits to the stellar continuum 
at red wavelengths and the excess continuum at bluer wavelengths. 



Based on similar work by IValenti et al.l (119931 ) and lHartigan et al.1 ()199ll ). we adopt initial 
values for our slab parameters of n=10 14 cm -3 , l=W 7 cm, and T=9000 K. In a few cases, the 
fitting routine had trouble establishing a good fit using our standard intial parameters. For these 
stars, the initial values were adjusted manually until a suitable location in parameter space was 
found. It should be noted that our criteria for a 'good fit' are mainly subjective, i.e. the fits to 
the data are examined manually. We look for two things: 1. Is the Balmer jump well matched 
by the slab?; and 2. Are the underlying photospheric features fit well by the chosen template? If 
these criteria are met and the model approaches the same reasonable reduced x 2 value for multiple 
models with different initial parameter values (this is usually an indication that the fitting routine 
has found a minimum x 2 i n the vicinity of the absolute minimum of the parameter space), the fit 
is considered 'good'. A discussion of the uniqueness of our models is given below. Typical reduced 
X 2 values for our models are ~1. 0-2.0. In order to produce superior fits, it was common to use a 
template that differed in spectral type by 1-2 subclasses from the determined spectral type of the 
star being modeled (see Table 2). This aspect of the model is discussed more thoroughly below. 

Initial values for f3 and Ay were determined by fitting the CTTS with a scaled template 
star in the absence of a slab. Due to the relative weakness of the excess spectrum compared to 
stellar photospheric emission at wavelengths > 4500 A (the exception being stars experiencing 
large amounts of mass accretion), this procedure generally established good fits to the red portion 
of the spectrum and provided stable starting points for the x 2 routine. Typical starting values 
for p were 10 -4 , corresponding to i?*~1.5 R for i?T~.75 R© and c?t~12 pc. Starting values for 
a were determined by assuming / aC c=-05. Again, some cases required manual adjustment of the 
initial values of the scale factors. 



The model slab flux is then generated with the best-fit parameters. The template photosphere 
is then subjected to the optical depth of the slab: 
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F T = [facce^ + (1 - facc)]F T . (B3) 

The slab is then scaled and added to the adjusted photosphere. The combined spectrum is then 
subjected to the reddening determined by the model and scaled to match the observed spectrum. 
The modeled flux can be summarized by the following equation: 

F* = (3 [aF slab + F^] (B4) 

where is the total (stellar plus accretion) flux, F s [ a t> is the slab flux, and F' T is the adjusted 
template star flux from eq. (B3). When constructing our models, we ignore any inclination or 
geometric effects concerning the location of the accretion column on the star. The slab is taken to 
be in front of the stellar photosphere along the line of sight to the observer. As a result, the filling 
factor we determine is the projected filling factor on the visible surface of the star. Ignoring these 
effects introduces some uncertainty in the final models. However, when averaged over the entire 
sample, assuming a uniform distribution of inclinations and number of accretion columns per star, 
this uncertainty should not affect the final correlations. 

We cannot ever be sure if our x 2 minimization routine approaches the absolute minimum x 2 
value for each fit. Instead, we can be confident in our final model parameters if, given a range of 
reasonable starting values, the fitting routine approaches a model of similar quality for each set of 
initial values. The uniqueness of the final stellar and accretion parameters is tested in the same 
manner. Figures Bl and B2 show this process graphically for the accreting star L1316. In these 
plots, all other intial parameters are fixed at the standard values described earlier in §4.1. Tables 
Bl and B2 list the initial parameters and final parameters for each model calculation. It is obvious 
that while the exact values of the slab parameters certainly vary with the values chosen for the 
intial parameters, the final stellar and accretion parameters do not change significantly for fits with 
similar x 2 values. Fits with the same x 2 value are practically indistintinguishable from one another. 
It is important to note that even though the final values for T and n vary by as much as a factor 
of 2 and as much as a factor of 10 for f3 and a factor of 5 for Ay, the optical depth through the 
slab, R*, and M do not vary appreciably at all. This shows that although the values of the slab 
parameters and scale factors are not unique to a specific x 2 value, the final stellar and accretion 
parameters are well constrained by the fitting procedure. 

We also investigated the effect of template spectral type on the final model parameters by 
fitting each spectrum with multiple template stars. The best fitting template is chosen manually 
by examining how well the underlying photospheric features are matched. It is sometimes the case 
that two different templates will yield similar stellar parameters and \ 2 values but one template 
does not appear to match the absoprtion features of the underlying star as well as the other. In 
these cases the template that visually shows the best match to the photospheric features is chosen 
for the model. More typically, however, changes in spectral type of the template by 1-2 subclasses 
result in a 10-30% change in the x 2 value for the fit. The exceptions to this procedure are L5924 and 
L4547. Due to our lack of a G-type template we modeled these stars using our earliest template, 
HD10476, a Kl star. Attempts to model these stars with later spectral types resulted in increasingly 



-29- 



worse fits. Surprisingly, fits to L5924 and L4547 using the Kl template produce decent fits without 
the inclusion of a slab. For this reason, as well as the fact that these stars do not show signs of 
significant mass accretion, they are grouped as NA stars and are not included in the disk-locking 
theory investigations discussed in §5. Thus, uncertainties due to the template spectral type, which 
for these two stars would be large, do not factor into our final analysis. 

The best-matching template for a star is usually quite obvious. The stars for which it is not 
so clear are those with higher accretion rates, and as a result more highly veiled spectra. However, 
these tend to be less sensitive to changes in template spectral type due to the dominance of the 
accretion emission. The stellar parameters for these stars tend to be more uncertain than for the 
stars with well-fit photospheres. Based on multiple fits to each star with different template spectra, 
we estimate that the model spectral types used are a ccurate to wi t hin 2 subtypes. From Table 2 
it is obvious that the spectral types determined by iRebull et al.l (J2002J) and the spectral types 
determined from our model fits do not a lways agree. I n fact , 30% of our best-fit spectral types 
differ by more than 2 subtypes from the IRebull et al.l (120021 ) spectral types. We do not i nclude 
L6175 in this number due to our lack of template beyond spectral type M3.5. IRebull et al.l (|2002j) 
estimate the uncertainties in their derived spectral types as 2 subtypes for K stars, and less than 
1 subcla ss for M star s . We attribute the differences in our model-determined spectral types and 
those of Re bull et al.l (I2002D to the fact that we are fitting a bluer portion of the stellar spectrum, 
whereas IRebull et al.1 (|2002j) used the region 5000 A< A < 9300 A to classify their stars. This 



spectral region is less affected by accretion emission and is particularly good for classifying M stars. 
We also do not have templates for each K and M subtype and so must use the closest available 
spectral type. Veiling due to accretion emission makes our spectra less sensitive to changes in 
template spectral type. We conclude that the differences observed here are not problematic to our 
analysis nor are they unexpected based on differences in the considered wavelength region. 
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Fig. 1. — Example spectra and their associated models. Panel (a) shows L1316, a typical SA 
star; panel (b) shows L3636, a typical WA star; panel (c) displays L5673, a typical NA star. 
The black lines are the observed spectrum; the final model spectrum is overplotted in red; the 
underlying template star spectrum is overplotted in green and the slab spectrum (with the Balmer 
lines excluded for clarity) is overplotted with a blue dashed line. Models are not produced for the 
NA stars. Spectra and model fits for the entire sample are given in Appendix A. 
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Fig. 2. — HR diagram of our sample using the lSiess et al.1 (|2000l ) PMS evolutionary models. The 
accreting sample is younger on average than the rest of the sample, as is to be expected if these 
stars are undergoing more significant mass accretion. The average age of our sample is ~6.4 Myrs, 
slightly older than the literature average for NGC 2264 of ~3 Myrs. Estimates of the uncertainties 
in T e ff and are shown above the legend. 



-8 -6 -4 -2 

log(L H/J /L ) 



Fig. 3. — The H/3 line luminosities of the SA and WA stars (black circles) plotted versus accretion 
luminosity. The best-fit line (solid) to the SA and WA stars is overplotted. The NA stars with 
measureable H/3 emission are plotted (red diamonds), though they are not included in the deter- 
mination of eq. (6). The NA star L6172 is plotted (open blue circle) using the values from |R02 
to determine an accretion luminosity. The vertical dashed line is a lower-bound on the amount of 
H/3 luminosity that can be attributed to accretion (see § 4.5). The bar at the bottom of the plot 
represents the range of H/3 line luminosity spanned by the WTTSs used to determine the average 
chromospheric contribution to the H/3 line flux . The equation of the best-fit line is almost identical 
to the relationship found bv lFang et al.l ([2009). 
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Fig. 4. — A comparison of Ay values determined from |R02| and those from this study. The line of 
equal values is overplotted with a solid line. There is poor agreement for ~60% of the stars. 
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Fig. 5. — A comparison of M values determined from IR02 and those from this study. The line of 
equal values is overplotted with a solid line. The |R02| accretion rates are systematically higher by 
an average factor of ~2.5. This factor tends to be higher for stars with lower accretion rates in our 
sample. 
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Fig. 6a. — Plot of eq. (1), predicted by the models of iKonigll (|199ll ) and IShu et al.l (j 19941 ). for the 
weakly accreting stars assuming a constant magnetic field from star to star. In Figures 6, 7, and 8 
the radius is in units of Rq, M* is in units of Mg, P ro t is in days, and M in units of 10 -7 M yr _1 . 
No significant correlation is present. 
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Fig. 6b. — Same as Fig. 6a but for the strongly accreting stars. The correlation is very weak. 
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Fig. 6c. — Same as figures 6a and 6b but for the combination of SA (blue circles) and WA (red 
diamonds) stars. Again, no significant correlation is found. 
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Fig. 7 . — Plot of eq. (1) for 26 stars from |R02| . The rotation periods are taken from iLamm et al 
(2004). There is no correlation present. 
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Fig. 8a. — Plot of eq. (3) for the weak accretors. The correlations for eq. (3) are much stronger 
than for eq. (1), with much smaller values of P. Overplotted is the line of best-fit (solid line) and 
the line of best-fit with the slope fixed at 1.0 (dashed line), as predicted by the models. 
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Fig. 8b. — Same as Fig. 8a but for the strongly accreti ng stars. The co r relat ion is strong with a 
very low value of P, suggesting support for the modified lOstriker fc Shul (j 19951 ) theory. 
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Fig. 8c. — Same as figures 8a and 8b but for the combined group of SA (blue circles) and WA (red 
diamonds) stars. This plot shows the strongest correlation with the best correlation significance. 
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Fig. 9. — A strong correlation is seen between M and B%, casting doubt on the significance of the 
correlation seen in Figure 8c. The associated correlation coefficient and correlation significance are 
given in Table 12. Plot symbols are the same as figs. 6 and 8. 
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Fig. 10. — Equation 11 plotted using the combination of the SA (blue circles) and WA (red 
diamonds) stars. A moderate correlation is found, although the correlation is increased significantly 
(r=.57, P=.004) if L7974, the lower right-most circle, is excluded from the analysis. The line of 
best-fit is overplotted. 
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Fig. 11. — Equation 12 plotted for the SA (blue circles) and WA (red diamonds) stars. No corre- 
lation is present, reaffirming the lack of support for the pure dipole models. 
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Fig. 12. — Theoretical equatorial dipolar magnetic field strengths at the stellar surface (upper 
panel) and rotation period (lower panel) plotted against the mass of the star's radiative core 



m 



( 


1999) for the 


Shu et al. 


( 


1994) 



constants included. Only the SA (blue circles) and WA stars (red diamonds) are included in the 
plot. Kolmogorov-Smirnov tests for the fully convective vs. partially convective stars yield high 
probabilities (Pxs , =-24, upper panel; Pks=-38, lower panel) that the samples are drawn from the 
same distribution. 
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Fig. Al. — Flux-calibrated spectra and model fits for our entire NGC 2264 sample. The SA stars 
are displayed first, followed by the WA and NA groups. The colors and line-types are the same as 
in Fig. 1. 
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Fig. Bl. — The effect of initial parameter choice on the final model values for L1316. The initial 
values of T and n are represented by green squares; the final values are shown with black asterisks; 
the values used in the final model are plotted with a blue upsidedown triangle. The dotted lines link 
the initial parameters with their final values. The final values lie along a straight line indicating 
the need for another constraint on the model. 
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Fig. B2. — Same as Figure Bl but for Ay and f3. Similar behavior is found for each pair of strongly 
coupled parameters. 
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Table 1. Observing log 



Observatory Night Conditions Slit width Instrument 



McDonald 28-Jan-2004 Photometric 4" LCS 
HJS 2.7m 29-Jan-2004 Light clouds 
30-Jan-2004 Photometric 



Kitt Peak 26-Feb-2005 Clear 3" RCS 

Mayall 4m 27-Feb-2005 Closed 

28-Feb-2005 Mostly clear 4" 

Ol-Mar-2005 Mostly clear 
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Table 2. NGC 2264 Sample 





Spectral Type 






ID 


Literature" 


Best-fit 6 


V-mag 


Group c 




M? 


M3.5 


"I 7 ^ 
1 / .00 


<5 A 


T A 1 8/1 


M2 


MO 


1 (\ R7 
10. D 1 




T A 1 Q9 


K4 


M2 


1 R 98 
10. Zo 






K6 


K7 


1 ^ ^/l 
10.0^1 


c: a 


T /1Q8fi 


Kl 


Kl 


1/11/1 
I't. 1*± 


<? A 


T Pil /I Q 


K4 


K7 


1 z. Q7 
10. o / 


<5 A 


T p^p^ 


M5 


M3.5 


1 Q 1 9 


^ a 


T PiP.7P. 


Ml 


K7 


i ^ no 
io. uy 


^ a 

DA 


L/ODoo 




K7 


i r or 

10. uo 


^ a 


T PiOOPi 

.L/oyuo 


K4 


K7 


1 ^ 07 
10. U / 


^ a 






K3 


1 ^ *}7 
10. o / 


<5 A 




M3 


M3.5 


1 7 ^ 

1 / .JO 


S A 






K7 




<5 A 


T £1 09 


K7 


M2 




c: a 


T £1 7Pi 


M7 


M3.5 


1 7 

i i .oy 


^ a 




K4 


K3 


10.00 


^ A 


T 1 70/1 


M2 


MO 


1 7 07 
1 1 .U / 


WA 


JjODoD 


M3 


M2 


1 7 RQ 

i ( .oy 


WA 

VV-TY 


-L/oDDD 


K7 


K4 


i r rr 

10.00 


WA 

VV-rV 




Kl 


K2 


1 R 0^! 


WA 
VV-rV 




K4 


K3 


1 ^ RO 

lO.OU 


AA/A 
VV-rV 


T /1^09 


K4 


MO 


io. oy 


WA 

VV-rV 


T /IQPi^ 


M3 


M2 




WA 


T Pil 08 


K7 


K7 


1 R 99 


WA 

VV-rV 


T £998 


Ml 


K7 


1 R 

io. oy 


WA 

VV-rV 


T /10Q8 


Ml. 5 


M2 


1 R 0/1 


i\r a 


L4443 


K4 


K3 


14.58 


NA 


L4511 




K7 


16.53 


NA 


L4547 


G3 


Kl 


13.57 


NA 


L5394 


Kl 


Kl 


14.70 


NA 


L5673 


K4 


K3 


15.39 


NA 


L5874 


Kl 


Kl 


14.36 


NA 


L5924 


G3 


Kl 


13.52 


NA 


L6024 




M2 


17.82 


NA 


L6172 


K5 


K3 


14.47 


NA 


L7714 


K4 


K2 


15.44 


NA 



a Spectral types determined by Rebull et al. (2002) 

b Spectral type of template star used to model the 
data. 

c SA=strongly accreting star, WA=weak accretor, 
NA=negligible accretor 
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Table 3. Main Sequence Templates 



Name Spectral Type B — V Rt Distance 









(R©) 


(PC) 


HD10476 


Kl 


0.84 


.86 


7.47 


HD109011 


K2 


0.94 


.83 


23.74 


HD45088 


K3 


0.97 


.82 


14.66 


GL570A 


K4 


1.11 


.79 


5.91 


GL394 


K7 


1.34 


.67 


10.99 


LTT11085 


MO 




.41 


30.48 


GJ393 a 


M2 


1.52 


.51 


7.23 


GJ273 a 


M4 


1.57 


.34 


3.80 



Note. — Rt is the radius of the template star. 
a Observed at Kitt Peak 



Table 4. Parameter Comparison 



K7 Main sequence template V830 Tau 



Star 


R, 


M 


face 


A v 




M 


face 


A v 




(Re) 


(10- 8 M yr- 1 ) 




(mags) 


(R©) 


(10- 8 M yr- 1 ) 




(mags) 


L4525 


1.59 


.87 


.021 


.05 


1.62 


.94 


.026 


.42 


L5108 


1.12 


.10 


.012 


~0 


1.21 


.10 


.006 


.46 


L5143 


1.70 


1.04 


.024 


.23 


1.66 


.85 


.023 


.48 


L5575 


2.03 


1.26 


.032 


.03 


2.11 


1.35 


.029 


.45 


L5638 


1.31 


.43 


.010 


.36 


1.12 


.32 


.038 


.52 


L5905 


3.44 


10.0 


.018 


1.08 


2.94 


5.56 


.028 


1.18 


L6039 


1.62 


1.42 


.019 


.32 


1.75 


1.76 


.017 


.62 


L6228 


1.22 


.22 


.014 


.93 


1.26 


.23 


.004 


1.16 



Note. — A comparison between the stellar and accretion parameters derived using a MS K7 
template and those derived using V830 Tau, a WTTS template. 
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Table 5. Stellar & Accretion Parameters - Strongly Accreting Stars 



Star 


M, 


K, 


L t 


Prot 


A v 


Age 


T eff 


face 


Lace 


M 




(M©) 


(Re) 


(Le) 


(days) 


(mags) 


(Myr) 


(K) 




(Lo) 


(10- 8 M yr- 1 ) 


L1316 


.16 


1.76 


.28 


4.55 


.21 


1.0 


3175 


.003 


.004 


.20 


L4184 


.40 


.88 


.14 


7.79 


.61 


7.0 


3800 


.028 


.069 


.73 


L4192 


.25 


2.54 


.78 


7.79 


1.33 


0.5 


3400 


.009 


.234 


11.4 


L4525 


.56 


1.59 


.58 


11.39 


.05 


2.0 


4000 


.021 


.064 


.87 


L4986 


1.85 


2.68 


4.36 


4.26 


1.18 


3.0 


5100 


.020 


.334 


2.32 


L5143 


.54 


1.70 


.66 


7.64 


.23 


3.0 


4000 


.024 


.069 


1.04 


L5355 


.12 


1.23 


.14 


4.17 


1.26 


1.0 


3175 


.003 


.016 


.79 


L5575 


.52 


2.26 


1.17 


11.73 


.16 


1.0 


4000 


.014 


.073 


1.26 


L5638 


.60 


1.31 


.40 


6.84 


.36 


4.5 


4000 


.010 


.041 


.43 


L5905 


.50 


3.44 


2.72 


8.46 


1.08 


0.4 


4000 


.018 


.301 


10.0 


L5967 


1.65 


3.15 


4.73 


10.77 


1.24 


1.0 


4800 


.089 


2.01 


18.4 


L6032 


.18 


3.57 


1.16 


.80 


1.31 


0.1 


3175 


.002 


.030 


2.83 


L6039 


.55 


1.62 


.60 


5.83 


.32 


3.0 


4000 


.019 


.101 


1.42 


L6102 


.22 


1.71 


.35 


9.04 


1.08 


1.5 


3400 


.012 


.023 


.87 


L6175 


.16 


2.02 


.37 


4.17 


.64 


0.5 


3175 


.001 


.008 


.48 


L7974 


1.60 


2.63 


3.30 


.88 


1.33 


2.0 


4800 


.099 


1.08 


8.53 
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Table 6. Stellar & Accretion Parameters- Weak Accretors 



Star 


M, 






Prot 


A v 


Age 


Teff 


face 


Lace 


M 




(M ) 


(Re) 


(Le) 


(days) 


(mags) 


(Myr) 


(K) 




(Lo) 


(10- 8 M yr- 1 ) 


L1704 


.40 


.72 


.10 


8.28 


.59 


11.0 


3800 


.037 


.005 


.04 


L3636 


.21 


1.33 


.21 


1.32 


.41 


2.0 


3400 


.001 


.002 


.06 


L3666 


.66 


.60 


.14 


11.21 


.01 


25.0 


4560 


.009 


.006 


.02 


L3748 


.92 


.99 


.53 


12.09 


1.53 


30.0 


4960 


.072 


.065 


.33 


L3809 


1.00 


1.18 


.66 


4.17 


1.30 


11.0 


4800 


.018 


.022 


.12 


L4602 


.40 


1.12 


.23 


6.97 


.01 


4.0 


3800 


.034 


.044 


.15 


L4956 


.21 


1.57 


.30 


5.51 


.49 


2.0 


3400 


.001 


.002 


.06 


L5108 


.61 


1.12 


.29 


6.40 


-0 


6.0 


4000 


.012 


.012 


.10 


L6228 


.60 


1.22 


.34 


8.28 


.93 


5.0 


4000 


.014 


.023 


.22 
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Table 7. Stellar Parameters-Negligible Accretors 



Star 


M t 


R, 


U 


A v 


Prot 


Age 


T eff 




(M ) 


(Re) 


(Le) 


(mags) 


(days) 


(Myr) 


(K) 


L4098 


.24 


1.35 


.35 


.32 


2.71 


.30 


3400 


L4443 


1.30 


1.32 


.94 


.27 


4.50 


5.0 


4800 


L4511 


.60 


1.06 


.26 


.01 


12.47 


7.0 


4000 


L4547 


1.30 


1.61 


1.57 


.06 


3.42 


8.0 


5100 


L5394 


1.20 


1.45 


1.28 


.64 


2.38 


10.0 


5100 


L5673 


.90 


.98 


.53 


.47 


8.28 


50.0 


4960 


L5874 


1.55 


2.05 


2.55 


.89 


3.70 


5.0 


5100 


L5924 


1.70 


2.31 


3.24 


.02 


3.00 


3.0 


5100 


L6024 


.22 


.80 


.12 


.56 


9.71 


1.0 


3400 


L6172 


1.60 


2.66 


3.37 


1.65 


1.76 


2.0 


4800 


L7714 


1.10 


1.30 


.80 


1.14 


4.45 


10.0 


4800 
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Table 8. WTTS H/3 line luminosities 



Star 


F Hf) a 

(10~ 14 erg cm -2 s _1 ) 


Lh/3 

(10- 5 L ) 


LkCa 3 


2.42 


1.48 


V826 Tau 


6.09 


3.72 


LkCa4 


5.86 


3.57 


V827 Tau 


6.22 


3.79 


V830 Tau 


4.08 


2.49 



Note. — H/3 line luminosities for the WTTSs 
used to determine the threshold in Figure 5 (dot- 
ted line). These stars show no evidence o f any 
near-IR excess (jKenvon &: Hartmannl 1 19951 s ) and 
only upper limits ar e found for the millimeter flux 
I Dutrev et al. 19961 ). The luminosities were calcu- 
lated assuming a distance of 140 pc to the Taurus 
star forming region. 

1 19931) . 



Valenti et al 



Table 9. Parameter Comparison 







R02" 








This study 




Star M. t 
(M ) 


log(Age) 
(yrs) 


R* Av 
(R Q ) (mags) 


log(M) 
(M yr" 1 ) 


M» 
(M ) 


log(Age) 
(yrs) 


(R Q ) (mags) 


log(M) Group 

(M© yr- 1 ) 



L4184 


.29 


5.6 


2.03 


.67 


-7.99 


L4525 


.68 


6.0 


1.93 




-8.24 


L4986 


1.46 


6.0 


2.66 


1.03 


-7.57 


L5143 


.65 


6.0 


1.90 


.61 


-7.72 


LI 704 


.35 


6.0 


1.56 


.43 


-8.73 


L3636 


.32 


6.2 


1.30 


.28 


-8.79 


L3748 


.74 


5.7 


2.43 




-8.25 


L3809 


.81 


6.4 


1.50 


.36 


-8.31 


L4602 


.65 


6.2 


1.56 


.44 




L4956 


.31 


6.0 


1.49 


.26 


-8.48 


L5108 


.66 


6.1 


1.71 


.98 


-8.61 


L4098 


.27 


5.2 


2.40 


.34 




L4547 


1.81 


6.7 


2.13 


1.03 




L5673 


1.21 


6.9 


1.38 


1.08 


-8.88 


L5874 


1.40 


6.4 


1.95 


.78 




L6172 


.91 


5.5 


3.21 


.57 


-6.89 



.40 


6.8 


.88 


.61 


-8.14 


SA 


.56 


6.3 


1.59 


.05 


-8.06 


SA 


1.85 


6.5 


2.68 


1.18 


-7.63 


SA 


.54 


6.5 


1.70 


.23 


-7.98 


SA 


.40 


7.0 


.72 


.59 


-9.37 


WA 


.21 


6.3 


1.33 


.41 


-9.21 


WA 


.92 


7.5 


.99 


1.53 


-8.48 


WA 


1.00 


7.0 


1.18 


1.30 


-8.90 


WA 


.40 


6.6 


1.12 


.01 


-8.81 


WA 


.21 


6.3 


1.57 


.49 


-9.24 


WA 


.61 


6.8 


1.22 


~0 


-8.98 


WA 


.24 


5.5 


2.58 


.32 




NA 


1.30 


6.9 


1.61 


.06 




NA 


.90 


7.7 


.98 


.47 




NA 


1.55 


6.7 


2.05 


.89 




NA 


1.60 


6.3 


2.66 


1.65 




NA 



Note. — A comparison between the stellar and accretion parameters derived by Rebull et al. (2002) and this study. Most 
mass measurements agree to within 30%. R02 calculates much younger ages for our sample. We attribute this to their use of 
the D'Antona & Mazzitelli (1994) stellar evolotion models which they show (R02 Fig. 13) predict ages much younger than 
those of Siess et al. (2000). 

a Rebull et al. (2002) 
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Table 10. 


Summary 


of the correlation 


results 




Eq. 1 




Eq. 3 






r P 


r 


P 


Slope 


WA stars 
SA stars 
All stars 


-.30 .43 
.26 .32 
.31 .13 


.88 
.91 
.91 


.002 

9xl0~ 7 
3xl0~ 10 


.57±.06 
.58±.06 
.63±.03 



Note. — Summary of the correlation results for the 
assumption of a constant magnetic field from star to 
star. Best-fit slopes for the poor correlations from eq. 
(1) are not computed. The slope of the best-fit line to 
the data plotted using eq. (3) is significantly different 
than the predicted slope of 1. r is the linear correlation 
coefficient; P is the correlation significance associated 
with the given value of r. 
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Table 11. Correlation results for M, i?*, and /< 









r 


P 


M 


vs. 


Rl 


.81 


8.2xlCT 7 


M 


vs. 


face 


.41 


.044 


M 


vs. 


JaccR* 


.80 


1.3xl0" 6 



Note. - r is the linear correla- 
tion coefficient; P is correlation sig- 
nificance associated with the given 
value of r. 
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Table Bl. Parameter tests for L1316: n vs. T 



Trial 


n 


To 


n 


T 


74000 




M 


x 2 




(10 14 cm- 3 ) 


(K) 


(10 14 cm- 3 ) 


(K) 




(R©) 


(10- 8 M yr- 1 ) 




1 


1.25 


9500 


9.47 


7150 


.21 


1.82 


3.12 


1.31 


2 


1.50 


8500 


8.56 


7470 


.21 


1.84 


3.15 


1.32 


3 


1.75 


10000 


7.97 


7455 


.21 


1.84 


3.13 


1.32 


4 


.75 


8000 


9.46 


7267 


.21 


1.83 


3.29 


1.31 


5 


5.00 


11000 


12.3 


6500 


.22 


1.75 


3.06 


1.30 


6 


3.00 


8000 


9.20 


7273 


.21 


1.83 


3.22 


1.31 


7 


3.50 


9000 


8.60 


7352 


.21 


1.82 


3.13 


1.31 


8 


4.00 


9500 


5.40 


8849 


.20 


1.88 


2.87 


1.34 


9 


4.5 


10000 


11.80 


6898 


.21 


1.80 


3.21 


1.31 


10 


1.00 


7000 


9.70 


7074 


.21 


1.81 


3.15 


1.31 


11 


1.50 


9500 


10.30 


7045 


.21 


1.81 


3.27 


1.31 


12 


1.75 


8000 


8.10 


7524 


.21 


1.85 


3.23 


1.32 



Note. — Initial and final parameter values for the density, n, and temperature, T. Each trial 
corresponds to a pair of points in Fig. Bl. 
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Table B2. Parameter tests for L1316: Ay vs. /3 



Trial 


Avo 


A) 


A v 




T4000 




M 


x 2 




(mags) 


(io- 20 


(mags) 


(io- 20 ) 




(Re) 


(10- 8 M yr- 1 ) 




1 


.15 


8.0 


.20 


6.65 


.21 


1.84 


3.13 


1.31 


2 


.20 


5.0 


.21 


6.75 


.20 


1.86 


3.18 


1.32 


3 


.25 


1.0 


.49 


9.08 


.36 


2.16 


5.94 


1.53 


4 


.30 


2.0 


.12 


6.00 


.22 


1.75 


3.19 


1.30 


5 


.40 


6.0 


.19 


6.54 


.21 


1.83 


3.11 


1.31 


6 


.50 


7.0 


.19 


6.54 


.21 


1.83 


3.12 


1.31 


7 


.70 


9.0 


.31 


7.57 


.19 


1.97 


3.76 


1.34 


8 


.85 


10.0 


.22 


7.11 


.27 


1.91 


3.25 


1.45 


9 


.90 


5.0 


.24 


6.51 


.17 


1.82 


4.08 


1.36 


10 


1.00 


4.0 


.49 


6.32 


.18 


1.80 


9.91 


1.87 


11 


.05 


9.0 


.07 


6.22 


.24 


1.78 


2.36 


1.40 


12 


.10 


1.0 


.10 


5.74 


.21 


1.71 


3.21 


1.33 


Note. 


— Same 


as Table Bl but for the flux scale factor /3 and the visual extinction 



Ay. These data correspond to Fig. B2. 



